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ABSTRACT 


This report details the development and provides the 
documentation for custom computer software applications that 
evaluate antenna parameters. The applications are written in 
Mathcad 3.1 for the following antenna types: linear, planar, 
and circular arrays; folded dipole; caged dipole; parabolic 
reflectors with helical and spiral feeds. Inputs to the 
applications are limited to the antenna’s physical dimensions. 
In some cases, ground parameters are required. 

The chapters are structured to provide the user with the 
necessary information needed to use and interpret the software 
for each antenna type. The software applications are provided 
aS appendices and give examples of each antenna type. 

Outputs of the applications provide various numerical and 
performance predictions, as well as far-field radiation 
patterns. The results computed are consistent with 
predictions provided in applicable publications, as well as 
those calculated by numerical antenna analysis programs such 


as ELNEC. 
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I. BACKGROUND AND PURPOSE 


This thesis details the development of user friendly 
mathematical applications capable of computing the radiation 
pattern and other pertinent antenna parameters of an antenna 
or antenna system based on available information, primarily 
dimensional information obtained from photographs. The 
applications developed are compatible with any IBM personnel 
computer using MS-DOS version 3.2 or higher and a math 
coprocessor. The applications are written in Mathcad version 
3.1 engineering software and can be run using either Mathcad 
3.1 or Mathcad version 4.0. These applications are intended 
for use by the Naval Maritime Intelligence Center 
(NAVMARINTCEN) . 

As already mentioned, user inputs are limited to antenna 
dimensions based on physical measurements primarily obtained 
from photographs, although other source data can be used. The 
applications are flexible to the extent that other parameters 
may be estimated and used as inputs to increase the accuracy 
of the computations to gain better insight into the antenna’s 
performance. The Mathcad applications provide various 
performance predictions as well as a graphical representation 
of the antenna’s far-field radiation pattern. The necessary 


background information needed to interpret the application’s 


formulas and displays are provided in the corresponding thesis 
Chapters. 

Dietrich [Ref. 1] and Gerry [Ref. 2] have completed the 
first and second reports of this project, respectively. This 
thesis is the last ina series of three reports intended to 


fulfill the NAVMARINTCEN statement of work. 


II. INTRODUCTION 


The increasingly sophisticated design by foreign 
countries of antennas warrants the development of improved 
analytical methods to provide timely and accurate technical 
assessment of these antennas. Through the process of reverse 
engineering, foreign communication, navigation, 
Identification, Friend or Foe (IFF), and radar antennas are 
evaluated to determine their capabilities, limitations, and 
vulnerabilities. Having limited information available 
concerning the parameters of an antenna, usually restricted to 
photographs that provide only the type and physical dimensions 
of the antenna, current methods of antenna evaluation are slow 
and tedious. 

With the availability of high-speed personal computers 
matched with sophisticated off-the-shelf engineering software, 
the assessment of antennas using only photographic information 
may now be obtained rapidly. The objective of this report is 
to provide NAVMARINTCEN with user friendly Mathcad software to 
achieve this type of performance in their evaluation of 
antenna systems. 

As with any computer analysis of antennas, the accuracy 
of the results depend upon the sophistication of software 
used, the complexity of calculations, the reliability of input 


data obtained from photographs, and the accuracy of any 


required input data that must be estimated. This report and 
the Mathcad applications are written with existing engineering 
equations for the type of antenna analyzed and inputs are 
limited to physical dimensions. Additional information such 
as ground effects and characteristic impedance may be included 
to increase the accuracy of the analysis. To this end, this 
report should provide an initial understanding as to the 
Capabilities of the antenna system in question. 

The chapters of this report document the equations and 
assumptions used for each specific antenna type and is written 
as a comprehensive reference for the software. Copies of each 
application are included as appendices to provide the user 


with a printed illustration of the software. 


IIt. THE UNIFORM LINEAR ARRAY ANTENNA 


Arrays are highly directional antenna system with narrow, 
steerable beams and low side lobes used for long distance 
communication and radar systems. The linear array is one 
geometrical configuration that consists of radiating elements 
lying along a straight line. As a result of the dependance of 
inter-element spacing (d) to the wavelength (A), linear arrays 
are inherently narrowband antenna systems. 

In this Mathcad application uniformly spaced, equally 
excited, isotropic point sources aligned on the z-axis are 
used as the radiating elements. In using isotropic point 
sources, polarization of the array is not calculated and an 
assumption is made that the antenna input resistance (R,,) is 
equa! to the radiation resistance (R_). This assumption 


results in calculations for gain (G) that are idealized 


values. The following inputs are needed to implement the 
software: 
N = number of isotropic radiating elements 
ad = inter-element spacing between adjacent elements 
6. = direction of main lobe “ 
f = frequency of interest “ 
I. = antenna feed current ~ 


N 
ll 


° e ° * 
characteristic feed impedance 


The first two inputs are obtained from visual data such as 


photographs, and the ~ indicates input parameters that are 
either known or are estimated. The linear array geometry is 


illustrated in Figure 3.1. 





FIGURE 3.1 Geometry of N-element linear array of 
isotropic point sources 


The total electric field of a linear array in the far- 
field, neglecting mutual coupling between adjacent elements, 
1s determined by the vector addition of the fields of 
individual radiating elements. It can be shown that the total 
field is equivalent to the product of the field of a single 
element, selected at a reference point, and a factor which is 


referred to as the array factor (AF): 


E(total) = [E(element)] x [AF] (aie) 


In (3.1), E is the vector electric field intensity, while AF 
is a scalar quantity. This result is referred to as pattern 
multiplication for arrays of identical elements. The array 
factor is a function of the geometry of the array and the 
element excitation amplitude and phase. For uniform 
excitation, by varying the inter-element spacing (da) and/or 
the phase (8) between the elements, the characteristics of the 
array factor, and consequently the total electric field, of 
the array can be controlled. [Ref. 3: pp. 204-207] 

This Mathcad application is written for an N-element 
linear array with equally spaced, uniformly excited isotropic 
radiating elements. In addition, element phase is assumed to 
vary linearly along the length of the array. The radiation 
pattern of a specific element type is neglected in (3.1) since 
in normal usage it will have little effect when the array 


consists of a large number of elements. However, the user of 


this software should be able to use information regarding the 
type of element and the orientation of the elements to narrow 
the range of possible 6.’s. The normalized array factor (AF) 


is fRef. 32: ppeewl2-2.4): 


sin( > (8) | 
AF,(0) = + (dimensionless) (3.2) 
a sin(> y (6) 
Z 
where 
w (8) = kdcos(68) + B (radians) ( 38 ) 
B = -kdcos(6,) (radians) ( 34 ) 


In (3.4), 9, is the direction of the maximum array Sealer 
value (1.e., for a broadside array 9, = 7/2, for an jeneeie 
array with the main beam directed at 180°, 0, = m7, and for a 
phased array with the main beam directed to +60°, 90, = 7/3), 
B is the phase excitation difference between the elements, k 
is the wavenumber (27/A), and dad is the inter-element spacing. 

With the linear array aligned along the z-axis, the array 
factor is a function of 8 but not of ¢. The visible region of 
the array in terms of 6, (8), and the inter-element spacing 


(d/A) is: 


0<O<n (radians) (33) 


B-kd < w(8) < Btkd (radians) (3.6) 


d/A determines how much of the array factor appears in the 
visible region as defined by (3.5) and (3.6). Although d/A 
provides guidance in determining the frequency to select for 
this Mathcad application, it should be noted that most arrays 
are designed with the inter-element spacing less than one 
wavelength, usually close to a half-wavelength. Thus, linear 
arrays are characteristic narrowband antenna systems. (Ref. 
ae. «6123 | 

An observation point is considered to be in the far-field 


if all of the following are satisfied [Ref. 4: pp. 24-25]: 


r21.6A (meters) (3.7) 
r 25 (Na) (meters) (3.8) 
je 2a (meters) (3%59)) 


The array length, Nd, in (3.8) and (3.9) is the maximum 
dimension of the antenna (the array length is assumed to 
include a distance of d/2 beyond each end element ([Ref. 5: 
p.55}]). The minimum distance to the far-field is found by 
selecting the largest value of (3.7), (3.8), and (3.9). 
Radiation intensity (U(8)) is related to the power 


radiated in a given direction and is independent of the 


distance to the observation point (Ref. 3: p. 27; Ref. 4: p. 
33]. For linear arrays U(8) is given by [Ref. 3: pp. 229- 
233): 


U(6) = (AF,)? (W/ solid angle) (3.10) 


Radiated power (P..,) is the total power radiated by the 
antenna and is obtained by integrating U(8) over a surface 
surrounding the antenna [Ref. 4: p. 33]. Since the geometry 
of the linear array is aligned on the z-axis and the array. 


factor is independent of ¢, P.., is thus [Ref. 3: pp. 220520ei% 


Pag = 27 [°u(8) sin(®) de (W) (3.1) 


Directivity (D,) is the maximum value of directive gain 
where directive gain is the ratio of U(8) in a specific 
direction to the average radiation intensity [Ref. 4: pp. 34- 
36). Hence, Dis given by? [Reta 3- pp ecsec one 


4n U. 


= ——__max (dimensionless) ( 3ez 


In (3.12), U,,, equals unity and occurs at 9, for the uniform 
linear array. 

Effective isotropic radiated power (EIRP) is a commonly 
used communications term that is defined as the product of 
antenna gain and total power accepted by the antenna from the 


transmitter. EIRP is determined as [Ref. 4: p. 62]: 
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D,  (W) (rein ) 


Assuming that the power delivered to the array consisting 
of N-isotropic point sources is equivalent to the radiated 
Peover, we Can estimate the radiation resistance (R.) of the 


antenna as [Ref. 3: p. 55]: 


peered eC) (3.14) 
[same 


Ir 


Using the same logic, we obtain the input resistance (R,,) of 
the array made up of isotropic point sources as equivalent to 
the radiation resistance: 


(Q) (iS 


The exact gain (G) of the array cannot be calculated with 
the data assumed, but an idealized value can be computed based 
on the assumption made to determine R,, and the assumptions 
made to calculate total antenna efficiency (€,). The total 


antenna efficiency is given by [Ref. 3: pp. 44-45]: 


€.=€,€.,=1-|C|? (dimensionless) (3-16) 

where 
Earl (dimensionless) (3.17) 
e =1-|I'|? (dimensionless) (3.18) 
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1 ee (dimensionless) (32259) 

Rin 7 Zo 
In (3.17), €,, 18 the conduction and dielectric efticremey 
which is unity for an ideal lossless antenna. In (3.19), @, 


is the characteristic impedance of the transmission line and 
Tr is voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(3.18) and the ideal gain for a linear array is [Ref. 3: pp. 
43-44): 


(dimensionless) (3.20) 


G( dB) =0 logy Gesn) (QB) (32.215 


As a result of constructing the linear array with 
isotropic point sources, the polarization loss factor (PLF) 
and antenna polarization are not calculated. Therefore, the 
maximum effective aperture (A,,) of the array when used as a 
receiver is calculated assuming negligible polarization 


mismatches. The A, is then [Ref. 3: p. 63]: 


A = GOA’ pre (m2) (3.22) 


em An 


where 


PLF =1 (dimensionless) (35223) 
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The assumptions used in finding R, and A, result in an 
feeaeuzed Value in Calculating the effective height (h,). The 


effective height is [Ref. 6: p. 42]):] 


faa | eee oe (meters) (3.24) 


eee) is the characteristic impedance of free space (1207). 

Two displays are produced in this Mathcad application. 
The first display is the rectangular representation of the 
magnitude of the array factor from -m to 7. The second 
display is the polar representation of the magnitude of the 
aitay Lactor. 

The Hansen-Woodyard end-fire array was not addressed in 
this application nor were arrays of uniform spacing but 
nonuniform amplitude (i.e., Dolph-Tschebyscheff array). In 
addition, nonuniformly spaced, uniformly excited arrays were 
not addressed by this application. 

The results obtained with this application are identical 
to those found in [Ref. 3: pp. 216-240} with respect to array 
factor patterns and directivity. Table 3.1 iS a comparison of 
Gata from (Ref. 3} for a 10-element equally excited, uniformly 
Spaced, end-fire linear array while Figure 3.2 is a comparison 
of the array factor pattern of the same array. Of note, the 
array factor patterns from this application and the patterns 


from {[Ref. 3] are identical as expected. 
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Table 3.1 End-fire Array Data Comparison 


ANTENNA REFERENCE CALCULATED 


DIMENSIONS DIRECTIVITY (D,) DIRECTIVITY (D,) 
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FIGURE 3.2 Array factor pattern of 10-element equally 
excited, uniformly spaced end-fire linear array with d = 


A/4m, A. = 0°, and B = -kd 


0 


> 


Table 3.2 is a comparison of data from [Ref. 3] for a 10- 
element equally excited, uniformly spaced, broadside linear 
array, and Figure 3.3 is a comparison of the array factor 
pattern of the same array. The array factor patterns from 


this application and the patterns from (Ref. 3] are identical. 


Table 3.2 Broadside Array Data Comparison 


ANTENNA REFERENCE CALCULATED 


DIMENSIONS DIRECrIVETY (D.) DIRECTIVITY (D,) 





16 





i 

5 

f 
Boye 


we 
\ 
\ 





FIGURE 3.3 Array factor pattern of 10-element equally 
excited, uniformly spaced broadside linear array with d = 
A/4m, 8. = 7/2, and B = 0 


Oo 


yy 


IV. THE PLANAR ARRAY ANTENNA 


Planar arrays, frequently used long distance 
communications and radar systems, exhibit characteristics 
analogous to those of linear arrays but have additional 
capability to control and shape the radiation pattern. Planar 
arrays are more versatile than linear arrays, providing 
symmetrical patterns with lower side lobes, and can be used to 
scan the main beam of the antenna array toward any point in 
space [Ref. 3: p. 261]. As with linear arrays, planar arrays 
are narrowband antenna systems due to the dependance of the 
inter-element spacing (d, and d,) on the wavelength (A). The 
following inputs are required analyze a planar array: 

M = number of isotropic radiating elements in the 

x-direction 

N = number of isotropic radiating elements in the 

y-direction 
d. = inter-element spacing between adjacent elements in 
the x-direction 
d, = inter-element spacing between adjacent elements in 
the y-direction 
= frequency of interest * 


direction of main lobe at 8 = 6.” 


f 
0, 


oO 


2, direction of main lobe at ¢ = @¢ 


I 


; antenna feed current “* 
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Z, = characteristic feed impedance ~ 


The first four inputs are data obtained primarily from visual 
data such as photographs, and the ~ indicates input parameters 


that are either known or estimated. The planar array geometry 


is illustrated in Figure 4.1. 





FIGURE 4.1 Geometry of M x N element planar array 
of isotropic point sources 
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This Mathcad application is written for an M x N planar 
array with uniformly spaced, equally excited, identical 
radiating elements positioned on the x-y plane. This 
application may also be used to approximate the array factor 
and field patterns for arrays that have elliptical or circular 
geometries by assuming an ellipse is equivalent to a rectangle 
or a circle is equivalent to a square, respectively [Ref. 6: 
p.187]. Since isotropic point sources are assumed as array 
elements, the polarization of the array is not calculated, and 
an assumption is made that the antenna input resistance (R,) 
is equivalent to the radiation” resistance Sikes This 
assumption results in values for gain (G) that are idealized. 
As with the linear array, the radiation pattern of a specific 
element type is neglected in the analysis of the planar array 
Since in normal usage it will have little effect for a planar 
array with a large number of elements. However, as with the 
linear array, the user of this software should be able to use 
information regarding the type of element and the orientation 
of the elements to narrow the range of the possible 0.,s and 
?,’'S- 

The normalized array factor (AF) is calculated for far- 
field observations neglecting mutual coupling between adjacent 
elements with equal amplitude excitation and is [{Ref. 3: pp. 


260-263]: 
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——_— > one 


sin( > y,(8,6) | sin( > v, (8,9) | 
A?) (4.1) 


mM six| oe | N sin| (8,4) 


where 
y,(68,6) = kd,sin(8)cos(@) + B, (radians) (4x2) 
ono kKarsin(e)sinth) + 6° (radians) (423) 
B, = -kd,sin(6,) cos(@,) (radians) (4.4) 
B, = -kKd,sin(6,) sin(@,) (radians) (Aiee5 ) 


Equation (4.1) 1s the normalized array factor that is a 
function of both 8 and @ where (4.2) and (4.3) are the array 
factor phase shifts in the x and y directions, respectively. 
(4.4) and (4.5) are the progressive phase shift between 
adjacent elements in the x and y directions, respectively, 


with the main beam directed along 9 = 0, and ¢ = @¢ Tora dees) 


& 
through (4.5), k is the wavenumber (27/A). 

When the inter-element spacing between the elements (d, 
and d,) equals or is greater than the wavelength (A), multiple 
maxima of equal magnitude are formed. To avoid creating 
multiple maxima, the inter-element spacing must be less than 


a wavelength at the frequency of operation [Ref. 3: p. 263]. 


With the dependance of inter-element spacing on _ the 


Zi 


wavelength, planar arrays are characteristic narrowband 
antenna systems. 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 


r2>1.6A (meters) (4.6) 


bo 5y (Ma,)* Steere (meters) (4.7) 


2((Ma)* + (Nay 


ry (meters) (4.8) 


Gers 
The length of an array is assumed to include a distance d/2 
beyond each end element [Ref. 5: p. 55). Therefore, the 


maximum dimension of the planar array is: 


j (Md,)? + (Nd,)* (meters) (4.9) 


The minimum distance to the far-field is found selecting the 
largest value of (4.6), (4.7), and (4.8). 
Directivity (D,) for planar arrays is determined 


approximately by [Ref. 7: p. 78]: 


AF_AF- 
D.=4n pAFn| 


O 


T= Pe aoe 4.10 
i i “(AF AF,)sin(@) a8 db ( ) 


The maximum value of the numerator in (4.10) is unity at 8 = 
6, and @ = ¢, since the array factor is normalized. The 


.e) 


denominator is the radiated power (P.,,) of the antenna where 
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P is obtained by integrating over all angles around the 


rad 
antenna above the x-y plane (0 < 8 < 7/2) [{Ref. 3: p. 28; Ref. 
7: p. 78]. This assumes the array is always pointed toward 
one side or the other of the x-y plane. 

Effective isotropic radiated power (EIRP) for planar 
arrays 1s defined as the preduct of antenna gain and the total 


power radiated. EIRP is determined as [{Ref. 4: p. 62]: 


EIRP = P,44D,  (W) (4.11) 


Assuming that the power delivered to planar arrays 
constructed of M x N isotropic radiating elements is 
equivalent to the radiated power, we estimate the radiation 


mearscance (R_) of the array as [Ref. 3: p. 55]: 


) 
a) (Api) 
eas 
Since the radiating elements are assumed to be ideal point 
sources, the input resistance of a planar array (R,,) is 
assumed to be equivalent to the radiation resistance: 


(QQ) (4.13) 


The gain (G) of a planar array cannot be calculated given 
the information assumed to be available, but an idealized 
value can be derived. Another factor affecting the gain is 
the assumption made when determining the total antenna 


efficiency (e,). The total antenna efficiency for an antenna 
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is (Ref. 3: pp. 44-45] 


6, ="bpeo, = 1 = |e (dimensionless) (4.14) 
where 
ej 2 (dimensionless) C41S) 
e.=1-|C|? (dimensionless) (4.16) 
ae ae . : 
Jhe See (dimensionless) (43470) 
I 5 wd 


In (4.14), €,, 1s the conduction and dielectric efficiency, 
which is unity for an ideal lossless antenna. In (4e7 jee 
is the characteristic impedance of the transmission line and 
ris the voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(4.15) when ideal lossless antenna is assumed. The ideal gain 
of a planar array in dimensionless and decibel quantities is 


then [Ref. 3: pp. 43-44]: 


G = €,23 (dimensionless) (4.18) 

G(aB) = 10 leq Ge 2) (dB) (4.195 

The maximum effective aperture (A,) for a planar array 

is derived by assuming polarization mismatches are negligible 


(this application is constructed using isotropic radiating 
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Oe See 


elements for which polarization loss factor (PLF) and antenna 


polarization are not calculated) [Ref. 3: p. 63): 


A - 20)" pir (m2) (4.20) 


em 4 


PLF = 1 (dimensionless) (4.21) 


The effective height (h,,) of the planar array is now 


estimated by [Ref. 6: p. 42}: 


Pe een C4222} 





In (4.22), n, 1s the characteristic impedance of free space 
(1207). 

Four two-dimensional array pattern displays are produced 
in this Mathcad application: the first is in the x-z plane; 
the second is in the y-z plane; the third display generated is 
a plane perpendicular to the x-y plane when @ = 7/4; the last 
display is a plane perpendicular to the x-y plane when @ = @4,. 
The last display is redundant when @¢ is either 0, 7/4, or 7/2. 

Results obtained with this application are the same as 
miaese found in [{Ref. 3: pp. 260-274]. Directivity was 
calculated using the approximate equation (4.10). An 
alternative approximation for directivity is given by either 


WRef. 3: p. 272-273]: 
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D, = Ween, len (dimensionless) (4223) 
or 
D, a (dimensionless) (4.24) 
QO, (rads?) 
In (4.23), D, and D, are the directivities of broadside @miie a. 
arrays of length and number of elements Md,, M, and Nd,, N in 
the x and y directions, respectively. In (4.24), fy, is the 
beam solid angle. 
Table 4.1 is a comparison of directivity computed for a 


10x10 planar array. The reference directivity is derived 


using (4.24) while the calculated directivity is from (4.10). 


Table 4.1 Planar Array Data Comparison 


ANTENNA REFERENCE CALCULATED 


DIMENSIONS DERECTINV ITY Sere) DIRECTIVITY “(ee 


23767 dae 24.07 dB 





Figure 4.2 and Figure 4.3 are comparisons of array factor 


patterns for a 5x5 planar array with patterns from [Ret 733) 
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and the calculated patterns. Figure 4.2 1s a representation 
of the array factor patterns in the x-z and y-z planes, and 
the scale of this figure has been increased in order to 
increase the clarity of the pattern. Figure 4.3 is the array 
pattern in the plane at @ = 45°, where the maximum beam 1s 
directed at ¢, = 45°. The array patterns produced by this 


application are identical to the patterns given in [Ref. 3]. 
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FIGURE 4.2 Array factor pattern of a 5x5 element, equally 
excited, planar array with d, = dad, = A/2 m, 8, = 30°,¢, = 


y oO 


45°, in the x-z and y-z planes 
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FIGURE 4.3 Array factor pattern of a 5x5 element, 
excited, planar array with da, =d 


45°, in the plane @ = 45° : 


Ix) 


equally 


= 4/2 m, 0, = 30°, @ = 


V. THE CIRCULAR ARRAY ANTENNA 


Circular arrays, or ring arrays, are antennas in which 
the radiating elements are placed on a circle with no elements 
positioned inside the circle. Circular antenna arrays with 
elements positioned inside the circle are planar arrays, and 
these antennas are addressed in Chapter IV. Circular array 
antenna are used for radio direction finding, air and space 
navigation, underground propagation, radar, and sonar systems 
(Ref. 3: p. 274]. As with linear and planar arrays, circular 
arrays with uniformly spaced, equally excited elements are 
narrowband antenna systems due to the dependance of the number 
of elements (N) and the circumference of the circle to the 
wavelength (A) [Ref. 3: p. 277]. 


The following inputs are required to analyze a circular 


array: 
N = number of isotropic radiating elements 
a = radius of the circle 
f = frequency of interest * 
6, = direction of main lobe at 6 = 6, 
@, = direction of main lobe at @ = ?, 
I, = antenna feed current “ 
Z, = characteristic feed impedance z 


The first two inputs are data primarily obtained from visual 


data such as photographs, and the ~ indicates input parameters 
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that are either known or estimated. The planar array geometry 


is illlustrated in Figure 5.1. 





FIGURE 5.1 Geometry of N-element circular array of 
1sotropic point sources 
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This Mathcad application is written for N-element 
circular arrays with uniformly spaced, equally excited, 
isotropic radiating elements on the x-y plane. The 
polarization of the array is not calculated since isotropic 
point sources are used, and the antenna input resistance @(Ry) 
is assumed to be equivalent to the radiation resistance (Rie 
This assumption results in the computation of ideal gain (G). 
As with the linear and planar arrays, the radiation pattern of 
a specific element type is neglected in the analysis of the 
Circular array since in normal usage it will have little 
effect for a circular array with a large number of elements. 
However, aS with both the linear and planar arrays, the user 
of this application should be able to use information 
regarding the type of element and the orientation of the 
elements to narrow the range of possible 0,’s and 9,’s. 

The normalized array factor (AF) is calculated for far- 
field observations neglecting mutual coupling between adjacent 
elements with equal amplitude excitation and is [Ref. 3: pp. 


2/4=278; Ref. 43° 350-254 18 


N 
AF(8,) = = @] [xasin@cos(9-9,) + %al (dimensionless 
n=1 
where 
® = 2n(—), n=1,2,...WN (radians) (5.2) 
a, = -Kasin@,cos(o,- ®,) (radians) (5.3) 


8/2 








Equation (5.1) is the normalized array factor. Equation (5.2) 
is the angular position of the nth element on the x-y plane. 
Equation (5.3) is the phase excitation of the nth element with 
the position of the main beam directed at 68 = 6, and ¢ = ¢,, 
and k is the wavenumber (27/A) and a is the radius of the 
array. 

An observation point is considered to be in the far-field 


if all of the following are satisfied [Ref. 4: pp. 24-25]: 


n> 16% (meters) (5.4) 
me 5 (2a ) (meters) (5.5) 
re atzay (meters) (57:6) 


The minimum distance to the far-field is found by taking the 
Maximum of (5.4), (5.4), and (5.6), where the diameter of the 
Circle is taken to be the maximum dimension of the antenna. 

fo determine the directivity (D,) for a circular array, 
the approximate equations of radiation intensity (U(8,¢)) for 
arrays and radiated power (P_,,) are used. Radiation intensity 
is [Ref. 3: pp. 229-233]: 


U(6,o) = (AF)? (W/ solid angle) (5.7) 


The radiated power is [Ref. 3: pp. 28]: 


Je 


ra 


ote */277(6,) sin(@) d@d> (Ww) (5.8) 
0 0 
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In (5.8), Prag 1S computed by integrating over all angles above 
the x-y plane (0 < 9 < w/2) [Ref. 7: p. 78]. This assumes 
that the array has a pattern maxima above the x-y plane with 
no radiation below the x-y plane. The directivity for the 


Circular array is [Ref. 3: pp. 229-233]: 


47 U 
oa (dimensionless) Boe} 
Pad 
In (5.9), Una 1S unity and occurs at 8 = 0, and ¢ = ¢, since 


the array factor is normalized. 
Effective isotropic radiated power (EIRP) is determined 
by [{Ref. 4: p. 62]: 


[256d a ee Oe) (5.10) 


Assuming that the power delivered to a circular array of 
N isotropic, uniformly spaced radiating elements is equivalent 
to the radiated power, we can estimate then the radiation 


resistance (R,) of the array as [Ref. 3: p. 55}: 


RL = — (Q) (5.11) 


Since the radiating elements are assumed to be ideal point 
sources, the input resistance of a circular array (R,,) is 
estimated to be equivalent to the radiation resistance: 


R, =R (Q) (5.12) 
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Since the input resistance is estimated, the gain (G) of 
a circular array cannot be calculated, but an ideal value can 
be derived. As with linear and planar arrays, addressed in 
Chapter III and Chapter IV, respectively, the total antenna 


efficiency (e€,) is [Ref. 3: pp. 44-45}: 


€. = €,€., = 1-|T/? (dimensionless) (5.13) 
where 
fae = 1 (dimensionless) (45) 
e.=1-|C |? (dimensionless) (5.15) 
ag > Be : . 
r= ——— (dimensionless) (5 2G) 
Rint Zo 
In (5.14), €.4 1s the conduction and dielectric efficiency, 
which is unity for an ideal lossless antenna. In (5.16), Z, 


1s the characteristic impedance of a transmission line, and [I 
is the voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(5.15) when an ideal lossless antenna is assumed, and the gain 
of a circular array is [Ref. 3: pp. 43-44]: 


G = e,D, (dimensionless) (San ) 


GiG@ee)) = i@wog,, (&,.D,) (dB) ( 52mike)) 
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The maximum effective aperture (A,,) of a circular array 
is derived assuming polarization mismatches are negligible 
(this Mathcad application is constructed using isotropic 
radiating elements for which polarization loss factor (PLF) 


and antenna polarization are not calculated) [{Ref. 3: p. 63]: 


em = G(A)* ppp (ne) (5.19) 
4 


PLF=1 (dimensionless) (5.20) 


The effective height (h,,) of a circular array is 


estimated by [Ref. 6: p. 42): 


bh. = Rr Aen (52233) 





In (5.21), n, is the characteristic impedance of free space 
(1207). 

Four two-dimensional array pattern displays are produced 
in this Mathcad application: the first is the x-z plane; the 
second is the y~-z plane; the third is a plane perpendicular to 
the x-y plane when @¢ = 7/4; the last is a plane perpendicular 
to the x-y plane when @ = @¢,. The last display is redundant 
when @ 1s either 0, 7/4, or 7/2. 

The directivity of a uniform circular array approaches 
the number of elements (N) as the radius (a) of the array 


becomes very large [{Ref. 3: p. 277]. Table 5.1 provides the 
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Girectivity of a uniform circular array as a function of the 


radius of the array. 


Table §.1 Circular Array Directivity Data 


az 1.59 


Figures 5.2 and 5.3 provide comparisons of the array 
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radiation patterns computed for a 10-element, uniformly 
spaced, equally excited circular array with a radius of 1.59 
meters with those in [{Ref. 3]. Figure 5.2 is the radiation 
pattern in the x-z plane, and the radiation pattern in the y-z 
plane is shown in Figure 5.3. For these patterns, 0, = @ = 
oF. The array radiation patterns generated by this 


application are identical to the patterns given in [Ref. 3]. 
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“0.8 -0.6 





FIGURE $5.2 Array factor pattern of a 10-element, uniformly 
spaced, equally excited circular array with a = 1.59 mn, 


and 8, = ¢, = 0°, in the x-z plane 
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FIGURE 5.3 Array factor pattern of a 10-element, uniform 
spaced, equally excited circular array with a = 1.59 nm, 
and 8, = ¢, = 0°, in the y-z plane 
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VI. THE FOLDED DIPOLE ANTENNA 


The folded dipole is a center-fed wire antenna that 
performs in a manner similar to that of a linear, center-fed 
thin-wire antenna with the added advantage of a more easily 
matched input impedance [Ref. 8: p. 7-36]. The folded dipole 
antenna consists of two parallel dipoles connected at the ends 
to form a narrow wire loop. The folded dipole antenna is 


illustrated in Figure 6.1. 





FIGURE 6.1 Geometry of folded dipole 
antenna 
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The folded dipole can be assumed to behave as a thin-wire 
dipole when the following are satisfied: the length of the 
antenna iS much greater than the diameter of the feed 
conductor (L > 2a); the distance between the conductors 1s 
greater than the radius of the feed conductor (d > a); 
(2Na/A)* < 1; all conductors are assumed to be perfect [Ref. 
Seep. 172); 

This Mathcad application calculates the parameters of a 
folded dipole by assuming the antenna behaves as a thin-wire 
dipole with an equivalent radius (a,). The current maximum is 
normalized to one amp. The equivalent radius is [Ref. 8: p. 


7-38]: 


2 2 
2 = exp a‘in(a) + b*in(b) + 2abin(a) (meters) (6.1) 


(a+b)? 


a 


This application provides the characteristics of the antenna 
both in free space and positioned horizontally over the earth. 
Although the folded dipole positioned vertically over the 
earth occurs frequently in VHF/UHF communications, it 1s not 
addressed in this research. 

The following inputs are required to analysis a folded 
dipole: 

a = radius of feed conductor 

b = radius of second conductor 
dad = spacing between conductors 
L 


= length of antenna 
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oy 
Il 


height of antenna above ground plane 


¥e 


f = frequency of interest 

€.4 = conduction/dielectric efficiency of conductors ~ 
€. = relative dielectric constant of ground plane " 

o = conductivity of ground plane * 


Yy, = Aistance for far-field parameter calculations 

Oo, = incoming wave electric field unit vector for antenna 
in free space 

Oo, = unit polarization vector for antenna in free space’ 

oh, = incoming wave electric field unit vector for 
horizontally positioned antenna ~* 

6h, = antenna polarization direction for horizontally 
positioned antenna * 

The first five inputs are physical dimensions obtained through 


photographs, the " indicates input parameters that are either 


known or estimated. 


A. FOLDED DIPOLE IN FREE SPACE 
Assuming the folded dipole behaves as a thin-wire dipole 
with a, < d < L, we have the far-field electric field 


intensity of the antenna aligned with the z-axis as [Ref. 3: 


peli2et. 
Te Jkt ee cos( > cos (8) = cos( == (6.2) 
JD Co) Sng), |) V/m ' 
vile PR Je sin (6) Ue 
In (6.2), n, 1s the intrinsic impedance of free space, equal 
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to 1207, and k = 2n/A is the wavenumber. The normalized 
electric field is found by dividing (6.2) by the maximum value 
of the electric field which occurs at 8 = 1/2 for L less than 
about 1.254. Clearly, the radiation pattern of the folded 
dipole is identical to that of a conventional dipole of length 
L. 

An observation point is considered to be in the far-field 


if all of the following are satisfied [Ref. 4: pp. 24-25]: 





> tsGa (mMebens) (673) 
ee Se (meters) (6.4) 
ae ai (meters) (6.5) 


The minimum distance to the far-field is found by comparing 
the values of (6.3), (6.4), and (6.5) and selecting the 
maximum value. The distance between the conductors is not 
taken into account in determining the minimum distance to the 
far-field since it is assumed that the folded dipole has an 
over all effective radius (a,). 

Directivity (D,) for the folded dipole is found by 
determining the radiation intensity (U(8)) and radiated power 
(P.a3)- Since the radiation pattern is symmetric in the x-y 
plane, U(68) and P,,, are functions of @ only. Radiation 


intensity is [Ref. 3: p. 28]: 
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Ps 
d= 5 |E(6)|*  (W/ solid angle) (6.6) 


The radiated power is [{Ref. 3: pp. 28]: 


Prag = 2n{"U(8) sin(®) de (W) (6.7) 
0 


Directivity for a folded dipole is [{Ref. 3: pp. 29]: 


4nU : 
i ins = (dimensionless) (6.8) 


rad 


where U,,, = U(7/2) . 
Effective isotropic radiated power (EIRP) is determined 
by [Ref. 4: p. 62]: 


D, (Ww) (6.9) 


The folded dipole is unique in its capability to act as 
a step up impedance transformer when the antenna is resonant. 
The antenna effectively operates as an unbalanced transmission 
line, and it can be analyzed by decomposing the current into 
two distinct modes: a transmission line mode [Figure 6.2] and 
an antenna mode [Figure 6.3]. fRef. 3: p. 341; Ref. 4: pp. 


205-207] 


44 








FIGURE 6.2 Transmission FIGURE 6.3 Antenna mode 
line mode 


The input impedance (Z,,) at the feed for the folded 
dipole is given in (6.10), and the radiation resistance (R,) 
for the antenna is the real component of (Z,,)) [Ref. 3: 


p.342): 


Ye EEE (fel) (6.10) 
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In (6.10), Z, is the input impedance of the transmission line 


mode and is obtained from the impedance transfer equation 


Z. = jz,tan{ <2) (Q) (6.15 


where Z, in (6.11) is the characteristic impedance of Sagtwas 
wire transmission line and is given in (6.12) [Ref. 3: pp. 


341-342; Ref. 8: p. 7-39]: 
7 =e 
Z, = —cosh Bea (Q) (6. Izy 
8 2 
The input impedance of a linear dipole of length L and 


equivalent radius a, is calculated by [Ref] 3: pp. ei2ayuere 


294, and 342]: 


7. = Rae ke aD (6.13) 
R 
Ry =) ee (Q) 
sin? aA Hooke) 
x 
Aa] = ; (2) 6.1 
sin a (6.15) 


: =! fy + In(kL) - C,(KL)] 


as) 
I 


No 
2 


+} 





[> sin( kL) (9,(2kD) - 28,(kL)]] (6.16) 


No 
ot 





1 kL 
5 cos (kL) (y ‘ 1n{ £2) : C,(2kL) -2.0;(kL)) 
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No 











X, = 7e[25;(kL)+cos(kL) [25; (kL) -S;(2kL)]] 
(en 175 
2 
BE sin( kz) {acj(kL)-¢,(2kL) ~cf 282 }) 
47 L 





In (6.16) and (6.17), C,(x) and S,(x) are the cosine and sine 


integrals, respectively, [Ref. 3: pp. 124 and 743-746]: 


Ge 


; = cos (t) = + — PA ely ae oe 6.18 
C, (x) if oer dt = y+ in(x) a ~ = ( ) 
Sax); = Sint) in(t) at = 2 1.) Sa oe Tani 6.19 


and ¥ 0.57721 (Euler’s constant). 

Mathcad computes summations more quickly than integrals. 
In this application, the summation forms of (6.18) and (6.19) 
are used with the index n ranging from 1 to 50 to ensure 
proper convergence. The maximum value of the argument of 
either C,(x) or S,(x) is x = 4mL/A, and for L = A the 
succeeding terms in the summation rapidly decrease for n > 20. 

The gain (G) of a folded dipole is obtained as the 
product of the antenna reflection efficiency (€,) and the 
directivity of (6.8) where the antenna reflection efficiency 


MompRef. 3: pp. 43-45): 


2, See, = pe) (dimensionless) (6.20) 
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where 


e.=1-|C |? (dimensionless) (6.21) 
Teele A e 

r= —— (dimensionless) (6222) 
Rint Zo 


In (6.21) e. is the mismatch efficiency and e is the 


r cd 


conduction and dielectric efficiency in (6.22). In (Ge 22) gem 
is the characteristic impedance of the transmission line given 
in (6.12). Therefore, the gain of the antenna is: 


G = e,D, (dimensionless) (6.23) 


G(dB) = 10log,,(€,D,) (dB) (6.24) 


The folded dipole in free space aligned on the z-axis is 
horizontally polarized with respect to the x-y plane. 
Therefore, given the direction of the incoming wave’s unit 
polarization vector (0,), the polarization loss factor (PLF) 


is [Ref. 3: pp. 48-53]: 


PLF = |G,° G;3|’ (dimensionless) (6.25) 


The maximum effective aperture (A,) of a folded dipole 


is (Ref. 3: p.63]: 


Aco St ee (m2) (6.26) 
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The effective height (h,,) of the folded dipole is (Ref. 


oop. 42); 


A Seen (mn) (6.27) 


The bandwidth of the folded dipole antenna is dependant 
on the length of the wires and has better bandwidth 
characteristics than a conventional dipole of the same length. 
It can be assumed that the bandwidth of the folded dipole is 
essentially the same as that of a conventional dipole of an 
radius a, (a < a, < ad/2) [Ref. 3: p. 346]. The bandwidth is 
increased when the distance between the conductors (d) is 
increased and/or when the conductor radii a and b are 
increased (Ref. 10: p. 190]. 

Tables 6.1 and 6.2 are comparisons of measured and 
calculated input impedance data for the folded dipole as a 
function of its electrical length. The radii of the 
conductors (a, b) are the same (0.0005 m) while the spacing 
between the conductors (d) is held constant at 0.00625 m and 


A = 1m [Ref. 8: pp. 7-40 - 7-41]. 
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TABLE 6.1 Folded Dipole Antenna Data Comparison 


ANTENNA MEASURED CALCULATED 
LENGTH (L) Re(Z,,) 2 Re(Z,) 2 
6533 ™ 


239 2 





100 
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TABLE 6.2 Folded Dipole Antenna Data Comparison 


ANTENNA MEASURED CALCULATED 


LENGTH (L) eZee Im(Z,;,) 2 





To further test the validity of the Mathcad application, 
a numerical antenna analysis program was also used. The 
program used was ELNEC, a powerful, easy-to-use program for 
modeling and analyzing a wide variety of antenna types 
including ground effects and parasitic structures. The 
fundamental computation portion of ELNEC is the same as 
MININEC (Version 3), that was developed by the Naval Ocean 
Systems Center [Ref. 11: pp. 5-7}. 

Table 6.3 is a comparison of the gain for a folded dipole 
in free space with different dimensions [Ref. 8: p. 3-29]. 


For each case, the radii of the conductors (0.0005 m) and the 


ayak 


distance between the conductors (12.5a m) are the same. As 
can be seen, the comparison between the gains obtained from 


various sources 1s very good. 


TABLE 6.3 Folded Dipole Antenna Data Comparison 


ANTENNA MEASURED ELNEC Mathcad 


DIMENSION Gain (G) Gain (G) Gain (G) 





The computed parameters in Tables 6.1, 6.2, and 6.3 are 
essentially the same as those measured or computed with ELNEC 
with exception of the ones marked with “". These discrepancies 
result in using the transmission line model and equivalent 
radius to analysis the folded dipole in this application. 

As previously mentioned, the folded dipole is a center- 
fed wire antenna that performs in a manner similar to a 
linear, center-fed, thin-wire antenna. Therefore, the 
radiation pattern of a folded dipole is also similar to that 
of a thin-wire antenna. The folded dipole radiation pattern 


in free space is illustrated in Figure 6.4. 
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FIGURE 6.4 Radiation Pattern folded dipole in free space 
with L = A/2 m, a = b = 0.0005 m, and a, = 0.004 m 


Figure 6.5 is a comparison of the radiation pattern of a 
folded dipole in free space obtained using ELNEC and the 
Mathcad application. The patterns are not identical, but are 
Similar and oriented in the same direction as expected. The 
difference can be attributed to the fact that ELNEC takes into 


account mutual coupling between the conductors while the 
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Mathcad application models the folded dipole as a thin-wire, 


center-fed dipole. 


=) 


~ Mothcod Pattern, ... ELNEC Pattern 





Figure 6.5 Comparison of radiation patterns using ELNEC 
and Mathcad of a Folded dipole in free space with L = A/2 
m, a = b = 0.0005 m, and a, = 0.004 m 
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B. FOLDED DIPOLE POSITIONED HORIZONTALLY OVER THE EARTH 
Through the use of image theory, the far-field electric 
field intensity of a folded dipole oriented horizontally above 
a flat earth can be obtained. The 8 component of the electric 
field of a horizontal dipole (parallel to the y-axis) in free 


space is (Ref. 3: pp. 143-144]: 


-jkryz| COS +2 sin (0) sing)| -cos( = 


Z | (V/m) 


6 
£0(8,o) =jn,—2— 


OUT ee ¥1-sin* (6) sin*() 


(6.28) 


From image theory, the 8 component of the electric field 


over ground is [Ref. 4: pp. 229-235]: 


E0,,, (8) = £0 (6, ) (e7khe0s ©) i Taam) (V/m) (6.29) 


In (6.28) and (6.29), E,(8,¢) is valid only above the ground 
(O < 6 < mw/2) (Ref. 3: pp. 135-142]. In (6.29), h is the 
distance from the earth to the antenna and the vertical 


Memlectionm Coefficient (f,) of the ground is (Ref. 4: pp. 231- 


Za3 }: 
/ _ — sine (0 
e’ cos(6) + e’,- sin? (6) 
where 
ge =€,- fj Tes NG dite Mad Olio cc) (6731) 


DS 


In (6.30) and (6.31), (€’,) is the relative complex effective 
dielectric constant of the ground, (e€,) is the relative 
dielectric constant of the ground, and (oc) is the conductivity 
of the ground. For a perfect ground plane, 0 > ~ and Yr, = l. 

The #@ component of the electric field of a horizontal 
dipole (parallel to the y-axis) in free space is [Ref. 3: pp. 


135-142]: 


) KL . kL 


EO G 70)! = 995 =e | (V/m) 


(6.32) 
2a 1-sin? (0) sin*() | 


Now from image theory, the # component of the electric field 


over earth is [Ref. 4: pp. 229-235]: 


Ed, (8) = EO), D(C aan aa a Pgem °o5) (V/m) (6.33) 


In (6.33), (T,) 1s the horizontal reflection coefficiene@ae- 


the ground and is [Ref. 4: pp. 231-232]: 


r. = cos(@) - /Je’,- sin? (6) 


SS ————————E (dimensionless) (6.34) 


cos (@) + /e’,- sin? (6) 


When CO. oO, ee = -l. 
The radiation intensity (U,(8,¢)) of a folded dipole over 
a earth must take into account both the 8 and @ components of 


the electric field. Therefore, the radiation intensity is 


Peet. 3: °p. 23am 
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2 
U,, (8) = 3 || 2Pcn (8) |? + |Zb,, (8) |?]} (W/ solid angle) (6-35) 





The radiated power (Ph,,,) follows (6.7) except that the field 


1s present only above the earth. Thus, 


Phyag = 20 {7 U,(8) sin(8) d® (Ww) (6.36) 


Directivity (Dh,) and effective isotropic radiated power 
(EIRPh) for a horizontal folded dipole over the earth are 
identical to (6.8) and (6.9), respectively, with the results 
of (6.36) inserted. 

Input resistance (Rh,,) for the folded dipole over the 
earth plane is found by (where I, in (6.32) is normalized to 


one amp) [Ref. 3: p. 124]: 


ee eee (6% 3m) 
a 


where the radiation resistance (Rh,) is, 


2 Ph 
(5, SEY) (6.38) 

2 

|Z 
The voltage reflection coefficient (IH), gain (Gh), 
maximum effective aperture (Ah,,), and maximum effective height 
(hh,,,) for a horizontal folded dipole over the earth are 


determined using the equations for those parameters for a 


folded dipole in free space, as listed above. 


a7 


The Polarization Loss Factor (PLFh) for the antenna over 
the earth can be determined when the direction of the main 
beam of the polarized antenna (8h,) is estimated. The antenna 
polarization vector (oh,) is then calculated, and estimating 
the unit polarization vector of the incoming wave (oh,), the 
PLF is found from (6.25). 

Table 6.4 1s a comparison of measured [Ref. 8: p. 3-29], 
ELNEC, and Mathcad software gains for a horizontal folded 
dipole over a reflecting sheet. The radii of the conductors 
(0.0067 m) are the same. As can be seen, the Mathcad software 


gain is similar to the measured and ELNEC gains. 


TABLE 6.4 Folded Dipole Antenna Data Comparison 


ANTENNA MEASURED ELNEC Mathcad 


DIMENSION GaumecG) Gain G) Gain sen 


A/2 m 


L/25.5 mi | 


A/8 m 
(a atl 





No input impedance data has been found to compare a 
folded dipole oriented horizontally over the earth with that 
computed by this Mathcad application. 

Figures 6.6 and 6.7 are the H-plane and E-plane radiation 


patterns of a folded dipole oriented horizontally over the 
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earth, respectively. Since the folded dipole is similar to a 
thin-wire center-fed dipole of length L, the illustrations 


below are similar to those of a dipole of length L. 
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FIGURE 6.6 H-plane radiation pattern of a folded dipole 
with L = A/2 m, a = b = 0.0005 m, dad = 0.00625 m, h = 0.5 
tie o = 0.01 S/m, e, = 15 
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FIGURE 6.7 E-plane radiation pattern of a folded dipole 
with L = A/2 m, a = b = 0.0005 m, da = 0.00625 m, h = 0.5 
m, 0 = 0.01 S/m, ¢, = 15 


Figures 6.8 and 6.9 are comparisons of the H-plane and E- 
plane radiation patterns between those obtained with ELNEC and 
Mathcad software of a folded dipole oriented horizontally over 
the earth, respectively. As expected, the patterns are 
Similar. As previosly mentioned, the difference can be 
attributed to the fact that ELNEC takes into account mutual 
coupling between the conductors while the Mathcad application 


is modeled after a thin-wire, center-fed dipole. 
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Figure 6.8 Comparison between ELNEC and Mathcad of H-plane 
radiation patterns of a folded dipole with L = A/2 m, a = 
oe O.0005 m, d = 0.00625 m, h=0.5 m, 0 = 0.01 S/m, te, = 
15 


Gz 0.4 


— Mothcod Pottern, ELNEC Pottern 





Figure 6.9 Comparison between ELNEC and Mathcad of E-plane 
radiation patterns of a folded dipole with L = A/2 m, a = 
be— 0.0005 m, ad = 0.00625 m, h = 0.5 m, 0 = 0.01 S/n, €, 
a5 
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VII. THE CAGED DIPOLE ANTENNA 


The caged dipole antenna is used in communications for 
its broadband antenna characteristics. Its performance is 
Similar to that of a thick cylindrical dipole. Since thick 
cylindrical dipoles are usually too clumsy to use except at 
very short wavelengths, thick dipoles are replaced with a 
caged dipole which consists of a number of thin wire dipoles 
arranged to form a circular cross section. The antenna is 
usually supported by a central conductor and support members 
that provide rigidity but contribute little to the electric 
properties of the caged dipole. The caged dipole (without the 
central conductor and the support members) is illustrated in 


Figure 7.1 (Ref. 12: pp. 172-174]. 


FIGURE 7.1 Geometry of caged dipole antenna 
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This Mathcad application models the caged dipole antenna 
as a circular array of center-fed, equally excited, uniformly 
spaced, thin-wire dipoles. This approach is taken to account 
for the number of conductors that are present in the caged 
dipole. As with the arrays previously considered, mutual 
coupling between the dipoles is neglected. In this manner, 
this application simulates the characteristics of a thick 
cylindrical dipole. The thin-wire dipoles are modeled after 
finite length dipoles in free space that have their length 
much greater than their diameter (L > wire diameter), and the 
dipoles are assumed to be perfect conductors. Since the 
central conductor and the support wires have little effect on 
the electrical properties of a caged dipole, mutual coupling 
from these components are not taken into account in this 
application. As seen in Figure 7.1, the inner ends of a caged 
dipole, the feed points, maybe coned in order to reduce base 
Gapacity [Ref. 12: pp. 172-174}. This effect is also not 
addressed. 

This application provides the characteristics of the 
antenna in free space as well as and oriented both vertically 


and horizontally over the earth. 


A. CAGED DIPOLE IN FREE SPACE 
The following inputs are required to analyze a caged 


dipole in free space: 
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N = number of conductors 
a = radius of antenna 


L = length of antenna 


f = frequency of interest ° 

I, = antenna feed current * 

Z, = input impedance * 

€., = conduction/dielectric efficiency of conductors * 
Y-, = distance for far-field parameter calculations “° 
Oo, = incoming wave electric field unit vector for 


ve 


antenna 


ve 


2s unit polarization vector for antenna 
The first three inputs are physical dimensions apeamea 
primarily from photographs, and the “ indicates input 
parameters that are either known or estimated. 

Assuming the caged dipole behaves as a circular array 
with equally excited (I,), uniformly spaced elements in the x- 
y plane with a radius of (a), we have the array factor 


(AF(6,@)) [Ref. 3: pp. 274-278]: 


N 
AF (0,0) = 3) aye” a ae (dimensionless) (7-1) 


n=1 


where the angular position in the nth element in the x-y plane 


is: 


®, = 2n(5), ee ee, ee (radians) (FZ) 
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Equation (7.1) is derived from: 


N 
AF (8,6) = >> Pe eee eC Gimensioniess) (7.3) 
n=1 
where the phase excitation of the nth element is: 


®, = 0 = -kasin6,cos(,-®,) (radians) (7.4) 


In (7.3), and the phase excitation of the nth element (a,) is 
zero as the result of desiring the main beam of the array to 
be directed perpendicular to the length of the antenna (8, = 
m/2). In (7.1), k is the wavenumber (27/A). 

Assuming each conductor of a caged dipole behaves as a 
thin-wire dipole, we have the far-field electric field 
intensity of a single conductor aligned with the z-axis as 


Mwtis 3s Pp. 120): 


Te ake cos(= cos (8)) - cos(=F] (7.5) 
E(0) =jn,————_ (V/m) ; 
2UI ep sin (6) 
In (7.5), n, = 120” and is the intrinsic impedance of free 


space. 

As previously discussed in Chapter III, the total 
electric field of an array constructed with identical elements 
is the product of the field of a single element and the array 
factor as a consequence of the principal of pattern 


multiplication PRet. 3... 207). Therefore, the total 
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electric field of a caged dipole in free space is approximated 
as: 


E,.(6,6) = £(6) AF(6,9) ( V/m) (7.6) 


An observation point is considered to be in the far-field 


if all of the following are satisfied [Ref. 4: pp. 24-25]: 


r21.6A (meters) (7an) 


r > Sv 2 eae (2a (meters) (7.8) 


2(L2 + (2a)?) 


Sa 1 


(meters) C7 2) 


The minimum distance to the far-field is found by taking the 
maximum of (7.7) through (7.9). In (7.8) and (7.9) ene 
maximum dimension of the antenna is the diagonal length of the 
array. 

Directivity (D,) for a caged dipole is found by 
determining the radiation intensity (U(@8,¢)) and radiated 


power (P,.,).- Radiation intensity 1s [Ref. 3: p. 28): 


2 
U(8,o) = aa IE.(0,@) |? (W/ solid angle) (Fino, 





The radiated power is [Ref. 3: pp. 28): 


2u fk : 
al ih U(8,@) sin(8) d6dpb (WwW) (7 mb) 
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BireceiVity for a caged dipole is (Ref. 3: pp. 29): 


4nU_., 


Bet (dimensionless) (12) 


rad 


where Un. = U(7/2,0). 


Effective isotropic radiated power (EIRP) is determined 


ere rt. 4:3: p. 62]: 


EIRP = P,,~D,  (W) CEeEe 


Input resistance (R,,) for a caged dipole in free space 


memmound from (Ref. 3: p. 124): 


R 





Roe eee Q 
in sin EE (£2) (7.14) 
2 
where the radiation resistance (R,) is: 
Pas & 
Ree SG) (Po 1US) 
jz," 


The gain (G) of a caged dipole is obtained as the product 
of the antenna reflection efficiency (e,) and the directivity, 


where the antenna reflection efficiency is (({Ref. 3: pp. 43- 


ars 


ye Ome C= €ey(1.- {TP |?) (dimensionless) (7 kG) 


where 


o. = ge) |e (dimensionless) (io dH) 
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| Oa a (dimensionless) (7.18) 


The parameter e, is the mismatch efficiency, e,, is the 
conduction and dielectric efficiency which is unity for a 
lossless conductor, and Z, is the characteristic impedance of 
a transmission line. Therefore, the gain of the antenna is: 


G=e,D, (dimensionless) (7.1393 


G (GB) =" 10 Nog.4 (ees (dB) (7220) 


A caged dipole in free space aligned parallel to the z- 
axis is horizontally polarized with respect to the x-y plane. 
Therefore, given the direction of the incoming wave’s unit 
polarization vector (o,), the polarization loss factor (PLF) 
is (Ref. 3: pp. 48-53]: 


—_ & 


PLF = |G, G3 |? (dimensionless) (7.21) 





The maximum effective aperture (A,,) of a caged dipole is 


[Rei . 3: peace. 


Zz 
Al. =] eee (m2) (7 a22) 
TU 


The effective height (h,,) of the caged dipole is [Ref. 


h = 2 fem (m) (7.29 














The bandwidth of a thick cylindrical dipole antenna is 
dependent on the ratio of the length of the antenna to the 
Giameter of the antenna (L/2a). When L/2a * 5000 the 
acceptable bandwidth of a thick cylindrical dipole is about 3% 
of the center frequency, and when L/2a *~ 260 the bandwidth is 
about 30%. Hence, if we assume a linear relationship between 
bandwidth and L/2a, and since the caged dipole’s performance 
1s Similar to that of a thick cylindrical dipole, then the 
bandwidth of a caged dipole is [Ref. 3: p. 333]: 

- f 


BW = f, (Hz) (7.24) 


igh low 


where the upper and lower frequencies are: 


eo a ee (Hz) (7.25) 
Pies ese Fisa}t (Hz) CDG) 
low Cc 22 Cc 


In (7.25) and (7.25), f, 1s the center frequency and F(L/2a) 
is the linear relationship function between antenna length and 
diameter. This relationship is obtained by solving y = mx + 
f. for m using x, = 5000, y, = 3% and x, = 260, y, = 30%. The 
relationship becomes: 


F{35)=-5.696-10° (35) +0-3248 (dimensionless) (7-27) 
a 
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Table 7.1 1s a comparison of the gain of a thick 
cylindrical dipole with that of a caged dipole (N = 8) in free 
Space computed both by ELNEC and Mathcad software [Ref. 8: p. 
2-38]. As previously mentioned in Chapter VII, ELNEC is an 
easy-to-use numerical analysis program for modeling and 
analyzing different antenna types [Ref. 11: pp. 5-7]. As can 
be seen, the comparison between the measured, ELNEC, and 


Mathcad application values are similar. 


TABLE 7.1 Caged Dipole Antenna Data Comparison 


ANTENNA MEASURED ELNEC Mathcad 


DIMENSIONS Gain (G) Gain (G) Gain (G) 





Figure 7.2 is the radiation pattern of a caged dipole (N 
= 8) in free space. It is essentially equivalent to that of 
a center-fed, thick cylindrical dipole having thegeeane 


conductor length and diameter [Ref. 8: p. 2-38]: 
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FIGURE 7.2 Radiation pattern of a caged dipole in free 
space with L = A/2 m, N = 8, and a = 0.125 m 


Figure 7.3 is a comparison of the radiation pattern of a 
caged dipole (N = 8) obtained using ELNEC and Mathcad. The 
patterns are not identical, but are similar and oriented in 
the same direction as expected. The difference can be 


attributed to the fact that ELNEC takes into account mutual 


vel 


coupling between elements while the Mathcad application 


ignores these effects. 
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FIGURE 7.3 Comparison of radiation patterns using ELNEC 
and Mathcad of a caged dipole in free space with L = A/2 
m, N = 8, and a = 0.125 m 
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B. CAGED DIPOLE ORIENTED VERTICALLY OVER EARTH 

The following known or estimated inputs are additional 
data required to analyze a caged dipole oriented vertically 
over the earth. The height of the antenna above ground is a 
physical dimension obtained through photographs: 


h = antenna height above ground 


€, = relative dielectric constant of earth 
Oo = conductivity of earth 
ov, = incoming wave electric field unit vector for 


vertical antenna 

8v, = antenna polarization direction for vertically 

orientated antenna 

As with the caged dipole in free space, the analysis of 
a caged dipole positioned vertically above the ground does not 
take into account the central conductor and support members 
Since they produce negligible effects on the antenna’s 
electrical properties. 

From image theory and pattern multiplication for arrays, 
the total far-field electric field intensity of a caged dipole 
oriented vertically above a flat earth is (Ref. 4: pp. 229- 
2a |: 


Ev,(6,6) =Ev(6) AF v (6, ) (e7kne0s 8) re Been) (V/m) C7 . 25) 


In (7.28), Ev(8) is the far+field electric field intensity of 
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a finite length dipole in free space which is given by (7.5), 
and AFv(68,¢) is the array factor for a circular array which is 
(7.1). The total electric field intensity is valid only above 
the ground (0 < 6 < 7/2) [Ref. 3: pp. 135-142]. The distance 
from the earth to the center of the antenna is (h), and the 
vertical reflection coefficient (T,) of the ground is [Ref. 4: 


Dp 2234s 2325r 


/ - Sa 
r,(6) = £7008 (0) See ae (dimensionless) (7.29) 


e/ cos (0) + fe’, -sin? (6) 


where 


(dimensionless) (72210 


In (7.29) and (7.30), (€’,) is the relative complex effective 
dielectric constant of the ground, (e,) is the relative 
dielectric constant of the ground, (oc) is the conductivity of 
the ground, and ¢, is the permittivity of free space (1/(367) 
x 10°). For a perfect ground plane, o + ~ and Tr, = 1. 

The radiation intensity (U,(8,¢)) of a caged dipole over 


earth is [Ref. 3: p. 28]: 


Z 





U,(8,o) = = (|Ev,(8,6) |] (W/ solid angle) (7.31) 


oO 


The radiated power (Pv,,,) 1S present only above the earth. 


Tus; [Re fe ass Sale 
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Pvraqg = { { ?Uv(8,6) sin(®) adh (Ww) Ciro) 
0 0 


Directivity (Dv,) and effective isotropic radiated power 
(EIRPv) for a vertically oriented caged dipole over the earth 
are given by (7.12) and (7.13), respectively, with the results 
Sieet7-52) inserted. Uvnax in (7.12) 1s simply the maximum 
value of (7.31) which, in Mathcad, is found by constructing 
(7.31) as a vector and then finding the maximum value of the 
vector. 

Input resistance (Rv,,), voltage reflection coefficient 
(TV), gain (Gh), maximum effective aperture (Av,,), and maximum 
effective height (hv,,) for a vertical caged dipole over the 
earth are determined using the equations for those parameters 
for a caged dipole in free space. 

The Polarization Loss Factor (PLFv) for the antenna over 
the earth can be determined when the direction of the main 
beam of the polarized antenna (@v,) is estimated. The antenna 
polarization vector (ov,) is then calculated and, estimating 
the unit polarization vector of the incoming wave (ov,), PLF 
memtOund from (7.21). 

Table 7.2 is a comparison of the gain of a vertical caged 
dipole over a flat earth obtained with both ELNEC and the 
Mathcad application. As can be seen, the gains are similar. 
As previously mentioned, the difference in the gains can be 


attributed to the fact that ELNEC takes into account mutual 
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coupling between the elements as well as ground effects. The 
Mathcad application ignores mutual coupling between elements 


and approximates ground effects. 


TABLE 7.2 Caged Dipole Antenna Gain Comparison 


ANTENNA ELNEC Mathcad 


DIMENSIONS GainmaG) Gainmaep 


0. 25. 


0.01 S/m 





Since the caged dipole behaves in a manner similar to 
that of a thin dipole except for bandwidth characteristics, 
the radiation pattern of a vertically oriented caged dipole 
Over ground will resemble that of a thin vertical dipole over 
ground. Figure 7.4 is a plot of the radiation pattern of a 
caged dipole over earth. There is a very close resemblance to 
the radiation pattern of a thin vertical dipole over ground 
(Ref. 4: p. 234]. Figure 7.5 is a comparison of radiation 
patterns computed by ELNEC and Mathcad software. As expected, 


the patterns are similar. 
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FIGURE 7.4 E-plane radiation pattern of a caged dipole 
with L = A/2 m, N= 8, a = 0.125 m, h = 0.25 m, o = 0.01 
S7m, and ¢, = 15 
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FIGURE 7.5 E-plane radiation pattern comparison between 
ELNEC and Mathcad of a caged dipole with L = A/2 m, N = 
Bea = 0.125 m, h = 0.25 m, o = 0.01 S/m, and e, = 15 
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C. CAGED DIPOLE ORIENTED HORIZONTALLY OVER EARTH 

The following known or estimated inputs are additional 
data required to analyze a caged dipole oriented horizontally 
over the earth: 

h = antenna height above ground 

€, = relative dielectric constant of earth 

o = conductivity of earth 

oh, = incoming wave electric field unit vector for 

vertical antenna 

6h, = antenna polarization direction for vertically 

orientated antenna 

As with the evaluations of a caged dipole both in free 
space and oriented vertically over the earth, the analysis of 
a caged dipole positioned horizontally above the ground does 
not take into account the central conductor and support 
members since they produce negligible effects on the antenna’s 
electrical properties. 

Through the use of image theory, the far-field electric 
field intensity of a caged dipole oriented horizontally above 
a flat earth can be obtained. The 8 component of the electric 
field of a horizontal dipole (parallel to the y-axis) in free 


Space is [Ref. 3: pp. 143-144]: 
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—— 


re okt cos(=-sin (8) sin@)} - cos 


EOh (8,0) =jn,————_ | (V 
. 2nT re v1-sin’ (6) sin*() on 


The 8 component of the total electric field over ground 


Mempneer. 4: pp. 229-235}: 


FON, (8,6) = £0h (8, ) AFh(8, d) (eskBeos (©) _ Te Skheos(®)) = (V/m) (7.34) 


The array factor (AFh(68,¢)) in (7.34) for a horizontal array 


is determined by [Ref. 3: pp. 274 and 776]: 


AFh(60,9) = y Fp se oe Sanat (dimensionless) (7-35) 
n=1 

In (7.33) through (7.35), the electric field (E6h(9,¢)), the 
total electric field (EO0h,(6,@)), and array factor (AFh(6,¢) ) 
are valid only above the ground (0 < 8 < 7/2) [Ref. 3: pp. 
135-142]. In (7.34), (h) 1s the distance from the earth to 
the center of the antenna, and the vertical reflection 

coefficient (I,) of the ground is given by (7.29). 
The ® component of the electric field of a horizontal 


dipole (parallel to the y-axis) in free space is [Ref. 3: pp. 


i142]: 
Te Jkt re cos(“sin (8) sin@)|-cos( | CTE) 
Eoh(8, 0) =jU,.—-——_— ——— a (V/im) ) 
eee 1-sin? (8) sin*(6) 


From image theory, the $ component of the total electric field 


over earth is [Ref. 4: pp. 229-235]: 
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Ebh, (6,6) =Bbh(8,) AFH(B, ob) (e7*Ros) + Pie Jkheos(9))  (y/m) (7.37) 


In (7.37), (AFh(8,¢)) is given by (7.35), and (f,) ise 
horizontal reflection coefficient of the ground and is (Ref. 
A: Pps [232i 


~ fe? —~ aine 
I’, _ £Co yee ee (dimensionless) 5 Be) 
cos (6) + ,/e’,- sin? (6) 


When o ~ ©, [, = -l. 

The radiation intensity (U,(9,¢)) of a caged dipole over 
a earth must take into account both the 8 and $ components of 
the electric field. Therefore, the radiation intensity is 
[Ref o35." pee. | 


Tee 
21. 





Uh (8,0) = [|Z0h, (8,6) |? + |Ebh, (8,0) \?] (W/ solid angle) (7.39) 
The radiated power (Ph,,,) is present only above the earth. 


Thus, 


Phyag = ff" Uh(8,6) sin(@) abd — (Ww) (7.40) 
Ira 0 0 


The remaining parameters (i.e., Dh,, EIRPh, Rh,,, etc.) 
for a horizontal caged dipole over the earth are calculated 
using the equations for these parameters listed previously. 

Table 7.3 is a comparison of the gain of a horizontal 


caged dipole over a flat earth obtained with ELNEC and the 
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Mathcad application. As can be seen, the gains are similar. 
As previously mentioned, the difference in gains can be 
attributed to the fact that ELNEC takes into account mutual 
coupling between the elements as well as ground effects. The 
Mathcad application ignores mutual coupling between elements 


and approximates ground effects. 


TABLE 7.3 Caged Dipole Antenna Gain Comparison 


ANTENNA ELNEC Mathcad 


DIMENSIONS Gain (G) Gain (G) 





As previously mentioned with regard to the vertically 
oriented caged dipole over the earth, the caged dipole behaves 
in a manner similar to that of a thin cylindrical dipole. 
Therefore, the radiation patterns of a horizontally oriented 
caged dipole are similar to those of horizontal thin 
cylindrical dipole which in turn resembles that of a thin 


dipole. Figures 7.6 and 7.7 illustrate the H-plane and E- 
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plane patterns, respectively, of a caged dipole over the earth 
(Ref. 4: pp. 231-2371 The radiation patterns of the 
horizontally orientated caged dipole over a flat earth are 
very similar to those of a horizontal thin dipole over earth. 
Figures 7.8 and 7.9 compare the radiation patterns between 
those computed with ELNEC and those computed by Mathcad 


software. As expected, the patterns are similar. 





FIGURE 7.6 H-plane radiation pattern of a caged dipole 
with L = A/2 m, N = 8, a = 0.125 m, h = 1.0m, o = 0.01 
S/m, and ¢, = 15 
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FIGURE 7.7 E-plane radiation pattern of a caged dipole 
with L = A/2 m, N = 8, a = 0.125 m, h = 1.0m, o = 0.01 
S/m, and ¢, = 15 


0.4 


Pattern 





FIGURE 7.8 Comparison between ELNEC and Mathcad of H-plane 
radiation patterns of a caged dipole with L = ijA/2 m, N = 
eeeea — 0-125 m, Hh = 1.0 m, o = 0.01 S/m, and ¢, = 15 
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FIGURE 7.9 Comparison between ELNEC and Mathcad of E-plane 
radiation patterns of a caged dipole with L = A/2 m, N = 
8, a = 0.125 m h=1.0m, 0 = 0.01 S/m, and ¢_, — ae 


84 


VIII. PARABOLIC REFLECTORS 


Many communications and radar operations require very 
large values of antenna gain (G) that are difficult to achieve 
with a single device. One of the most popular ways to 
construct a very high gain antenna is to reflect the 
electromagnetic energy of a relatively small antenna off a 
metallic parabolic dish in the far-field of the feed device 
(Ref. 13: p. 12-2]. If the source antenna is located at the 
focus of the parabola and if the radius (a) of the parabolic 
dish periphery is large with respect to wavelength (A), then 
the reflector may be approximated as a uniform phase aperture 
antenna. Gain in excess of 30 dB is common for this type of 
structure [Ref. 4: pp. 423-424]. The geometry of a parabolic 
reflector is illustrated in Figure 8.1, where the parameters 
include the focal distance (f), the distance (p) from the 
origin to any point (P) on the reflector surface, and the 
distance (r’) from the origin to the projection of point (P) 
onto the z = O plane. The primed angles in Figure 8.1 
correspond to the parameters associated with the feed horn. 
The Mathcad application user should note that the focal 
distance should exceed the minimum distance to far-field for 


the feed antenna. 
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FIGURE 8.1 Parabolic Reflector Geometry 


The Mathcad parabolic reflector applications are unique 
in that the software is not written as a stand alone package. 
Rather, the reflector software is run as an addendum to the 
programs that already exist for various antennas [Ref. 2: pp. 
4-15, 67-114}. Four types of feeds are analyzed: helical, 


planar equiangular spiral, pyramidal horn, and conical horn. 
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The following additional inputs are required to analyze 


parabolic reflectors: 


foc = reflector focal length 

aes = radius of the mouth of the reflector 
E, = electric field scale factor " 

€, = blockage efficiency * 


Be — Spar efficiency ° 


Yero = Secondary far field observation distance “ 
X,yY,Z = polarization loss factor coordinates “* 
Oo, = incoming wave unit polarization vector * 
tli, t2 = radiated power increments for @ and 8 
il = secondary field increments 

N,M = summation increments 


The first two inputs are physical dimensions obtained 
primarily through photographs, inputs with the “ are 
parameters that are either known or estimated, and the last 
inputs are only used to affect the processing time of the 
application. 

In order to calculate the far-field patterns of a 
parabolic reflector, the electric field at the aperture of the 
device must be known. Since the electric field patterns of 
interest are in the far-field, it is assumed that the circular 
aperture created by the reflector is located at the origin and 
is parallel to the z = 0 plane. The Cartesian components of 


the aperture electric field are expressed as [Ref. 14: pp. 


mpea 2s; Ref. 15: pp. 29-31): 
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Eax(0',9’) = [- Ep, (0' (4) , 9") cos(¢’) ae 
x yet) (r’) , 0’) sin(¢’) ] eo eee (V/m) 


E,y(0/, 0’) = [Ege(8/(r’) 9’) sin(¢’) an 
+ Ey, (8! (r/) 16/) Cosegi me ies 272) amma ats 


In (8.1) and (8.2), (E,,(r ,8')) and (E,(r 9G) ) eee 
aperture electric fields and k is the wavenumber (27/A). 
The relationships between (89’), (p), and (r’) needed to 


evaluate (8.1) and (8.2) are expressed as [Ref. 4: p. 426]: 


i 
Q/ = 2tan'i{ =| (radians) (8.3) 


4 ( foc)*+(r’)* 


8.4 
4foc 0) ! 


p(zr’) = 

The electric field from the entire parabolic assembly is 
commonly referred to as the secondary pattern. Referring to 
Figure 8.1, the electric field in the -z half-plane is 
considered the secondary pattern. The secondary field of a 
parabolic reflector is calculated using the electric field 
integral equation solution and computing the component vector 


potentials (P,(68,¢)), (P,(8,¢)) as follows [Ref. 4: p. 426]: 


P(0; d) -{"f Aref Ee aa ’) [Geet ca cos (o) + Sind vena) 
“ae (8.5) 


ig ole wen (V/m) 
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rs) ye es / / jkr’sin(@)(cos(¢’ )cos(¢) + santa rainc99)| 
P5(8,9) =|0 [E(x 9’) le 
pO (8.6) 


xr/ dar’ ag’ (V/m) 


Double integrals in Mathcad take significantly longer to 
converge to an answer than do single integrals [Ref. 15: p. 
200}. To reduce the processing time of the applications, 
(8.5) and (8.6) were each separated into two single integrals. 
Therefore, for the applications, (8.5) and (8.6) are each 


computed as: 


P(0,) = [°"A(0,6,6/) do’ (V/m) (8.7) 


where 


A(6,, db’) = f° E. (r/, ’) ae jkr'sin (8) [cos ($’) cos() + sin(o’)sin(o)] -/gy/ (8.8) 


Any contribution from feed system back lobes is ignored, 
and the component vector potentials calculated in (8.7) and 
(8.8) are applied to the subsequent formulas to compute the 


secondary field [Ref. 4: p. 383]: 


: -IKI gz 
E, (8,0) - Jke "1p, (8, ) cos () + P,(8,o) sin() ] (V/m) (8.9) 
OUI reo 
: We 7 
E,(8,) = =———_—_[-P, (8, @) sin(o) + P,(8,) cos () jcos (8) (8.10) 
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The radiation intensity (U(8,¢)), radiated power (P,.,), 
and directivity (D,) of the secondary field pattern are 
expressed using the intrinsic impedance of free space (n, = 


1207) as follows [{Ref. 3: pp. 27-30]: 


2 
U(8,) = “Site |B (8,4) P+ [Ey (8, ) |7] (W/solid angle) (8.11) 
2n rn/2 ; 
eee =f I, U(8,) sin (8) dB dob (W) (8.12) 
4nU 
o = nex (dimensionless) (Sa 1a) 
Pons 
In (8.13), (U,,,) 1s the maximum value of radiation intensity, 


obtained by creating (U,,,) aS a matrix and finding the maximum 
value of the matrix. 

The computer processing time required to evaluate (8.12) 
is prohibitively long. Thus, the parabolic reflector Mathcad 
applications use the following numerical approximation to 


evaluate radiated power [Ref. 3: pp. 37-42]: 


N M 
_ {27 us 8.14 
Prad = oule” bs U(G,, O,) SimcGR) (W) ( ) 
In (8.14), (N and M) are user selected integers. This 


approximation for P,,, reduces the processing time of the 
applications substantially. As N and M increase, the results 


obtained with (8.14) converge to that obtained with (8.12). 
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When N = M = 3, the processing time for only P,,, of the 
helical feed reflector option is approximately four hours on 
a 33 MHz, 386 machine, and when N and M is increased to 5, the 
time to calculate P,,, for the same application is ten hours. 
For the applications, we chose N = M = 10 to evaluate P,,, with 
reasonable accuracy. With this choice, it takes approximately 
one day for the application to run, but this is much less than 
the thirty-six hours it takes if (8.12) is used. 

Parabolic reflectors are members of the aperture antenna 
family. The gain of an aperture antenna is the product of its 
feemeemne efficiency (ce,,) and directivity. Therefore, the gain 
of a parabolic reflector is written as [Ref. 6: p. 575): 


eG (dimensionless) (8.15) 


i 
fer] 
. 
o 


Aperture efficiency is the product of several separate 
terms, including spillover efficiency (e,), taper efficiency 
fee Phase (or random surface error) efficiency (e€,), 
polarization efficiency (€,), blockage efficiency (¢€,), spar 
Emimetency (¢.,), and’ohmic efficiency (€¢f;,) - 

Spillover efficiency (e,) is the fraction of power 
radiated by the feed system that is intercepted and reflected 
by the parabolic dish. Spillover efficiency is the most 
important efficiency term for any feed system and is defined 


See [Ref. 6: pp. 583]: 
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("rae sin (6) dd’ dd’ 
= (dimenstontess) (8.16) 

[" f°a,(8 sin (8) de’ ag’ 

0 0 


In (8.16), (U,;(9’)) is the radiation intensity function of the 
feed antenna and (68,) is the value of (8.3) when (r’) equals 
the radius of the mouth of the reflector (a,,;). 

Taper efficiency (€,) accounts for the lack of amplitude 
uniformity of the feed pattern on the surface of the 
reflector. Taper efficiency, which is a dimensionless 


quantity, is expressed as [Ref. 3: pp. 626-627}: 


Tr / / 
[ [his Ane can| © lc6' ag’ 
0 Jo Jager 2 


2n 6, 4nU, (8 : / / / 
f. i. pi nee 





c,22cot?{ Se) (8.17) 
Z 


Phase efficiency (e€,) 18 a function of the "powerwaiee 
that occurs if the field at the mouth of the dish is not in 
phase at every point in the aperture. Phase efficiency is 
approximately given by [Ref. 4: pp. 434-435]: 


-4n6°10°4,.¢ 


: (8.18) 


(dimensionless) 


The power that is lost due to cross-polarized electric 
fields in the aperture of the reflector determines 


polarization efficiency (e,). Polarization efficiency 1s 


SZ 


g@ifficult to precisely calculate, but can be reasonably 
approximated as 0.98 [Ref, 4: p. 435]. 

Blockage (€,) and spar (€,,) efficiencies are determined 
by the presence of feed structures, support struts, and 
possible signal processing equipment in front of the 
reflector’s aperture. Precise computation of blockage and 
spar efficiencies requires an extensive method of moments 
simulation and is not attempted by the Mathcad parabolic 
reflector applications. However, if the radius (a,;) of the 
blocking structure at the focal point and the number of spars 
(N) that are (A/2) thick is known, blockage efficiency and 
Spar efficiency can be estimated using Tables 10.1 and 10.2, 


respectively (Ref. 4: p.436]. 


TABLE 8.1 Blockage Efficiency 





3) 6) 


Ohmic efficiency (€ jm.) for a reflector is usually very 
high but difficult to accurately calculate. The Mathcad 
parabolic reflector applications assume ohmic efficiency is 
O23. 

With all reflector efficiency terms known or estimated, 
the aperture efficiency (¢€,,), effective aperture (A,), and 


isotropic radiated power (EIRP) are given by [Ref. 3: pp. 623- 


630i]; 
bans Eso e ope rene sole (dimensionless) (82 Eas 
Ag = feet iases) (m*) (8.20) 
EIRP = P_,gD, (WwW) (8.21) 


The magnitude of the feed current (I,) for the source 
antenna cannot be estimated based on measured geometry. 
However, if the feed current is normalized to 1 amp, radiation 
resistance (R,) for the entire reflector structure is: 


Tie 





rad 
= (Q) (8.22) 
r 
eis 
In (8.22), (P..44) 1S the radiated power from the aperture of 


the reflector, not of the feed antenna. 
The Mathcad parabolic reflector applications use the 


radiation resistance computed for the parabolic reflector 
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antenna system to determine effective height (h,,) by [Ref. 6: 


p. 42]: 


puto |e (Ses) 





The bandwidth of a parabolic reflector is assumed to be 
the same as its feed antenna. 

The polarization loss factor (PLF) of the parabolic 
reflector for a user selected point in the secondary far-field 
and a specific wave unit polarization vector (0,) is found 
using the electric field components of (8.9) and (8.10) as 


momemowsS (Ref. 3: 51-53]: 


PLF = |¢,-G,|? (dimensionless) (8.24) 


where the unit polarization vector of the antenna is: 


4,E,+ 4,E,+ 4,£, 


Vee yall 


Cmex, V,Z) = (dimensionless) (8.25) 


and the Cartesian components of the electric field are [Ref. 


Seem. | 76); 
E,= Eycos (8) cos() - E,sin() (V/m) (8.26) 
E,,= Eycos (8) sin() + Eycos (o) (V/m) (S27 ) 
E,= -E,sin(@) (V/m) (8.28) 


o> 


All of the secondary field parameters calculated by the 
Mathcad parabolic reflector applications assume that the 
observation point is in the far-field. Thus, all of the 
subsequent formulas must hold for the application calculations 


to be valid [Ref. 4: pp. 24-25}: 


D2 1 eon (m) (8.29) 
it 2 Care (m) (8.30) 
fe > —— (dimensionless) (8.3) 


Table 8.3 iS a comparison of reflector directivity 
computed by the Mathcad applications to predictions made for 
an identical antenna with unity aperture efficiency [Ref. 3: 
Dee “6547 As can be seen, the comparisons vary by 
approximately 8 dB. This can be attributed to the fact that 
P.1qg WAS approximated using (8.14) instead of (8.12), and also 
to the fact that the numerical tolerance variable (TOL in 
Mathcad’s math menu) for the application was increased to 0.5 
verses 0.001 in order to reduce processing time. This 
increase in TOL reduces both the accuracy of (8.14) and the 


processing time of the application. 
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TABLE 8.3 Parabolic Reflector Data Comparison 


a 


ANTENNA PREDICTED CALCULATED 


PARAMETER DATA DATA 


DIRECTIVITY 42 @B | 50.2 aB 





Figure 8.2 1S a comparison of measured data fed by an 
optimum horn with that calculated by the Mathcad application 
fed by an equiangular spiral for a parabolic reflector (a,,.; = 
0.320 m, £f = 0.147 m) operating at 35 GHz [Ref. 14: p. 126]. 
As can be seen, the comparison 1s good. The Mathcad pattern 
of Figure 8.2 is not symmetric which can be attributed to the 
number of increments used to increase the speed of processing. 

The reflectors fed by pyramidal and conical horn antennas 
are not included in this report. The applications have been 
developed and are included on the diskette containing the 
software. The complexity of the software results in 
prohibitively long processing time for the two applications 
(1.e., the pyramidal horn fed reflector application was still 
calculating P,,,, after seven days when run on a 50 MHz, 486 


personal computer using Mathcad 4.0). 


oF 


-0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 0.1 


~ Mothcod Pattern, ... Measured Pattern 





FIGURE 8.2 Comparison of Parabolic Reflector E-Plane 
Electric Field Patterns 
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IX. REMARKS AND CONCLUSIONS 


This thesis and the accompanying software applications 
are intended to provide the Naval Maritime Intelligence Center 
(NAVMARINTCEN) with software for the analysis of various 
antenna types where information regarding the antenna is 
restricted primarily to the antenna type and the physical 
dimensions of the antenna. The goal of this research was to 
provide user friendly software that an engineer familiar with 
basic antenna types can easily use and interpret. Although 
the software was developed for that purpose, the complexity of 
formulas that apply to many of the antenna types of necessity 
reduced the simplicity of the applications. 

The software applications developed are written in 
Mathcad 3.1 and are compatible with any IBM personal computer 
that supports Mathcad 3.1 or Mathcad 4.0 for Windows. Many of 
the applications require extensive processing time. In order 
to reduce processing time, which results from the complexity 
of various integrals, approximate numerical techniques were 
implemented to evaluate the integrals and the numerical 
tolerance of the more complicated applications was reduced. 
Even with the processing time reduced, some programs require 
up to five days to complete when run with Mathcad 4.0 ona 50 
MHz, 486 personal computer. Applications that did not exceed 


a day to complete were run on a 33 MHZ, 386 personal computer 


29 


with Mathcad 4.0. All applications were run using Mathcad 4.0 
Since Mathcad 4.0 1S approximately twice as fast as Mathcad 
3.1. AS previously mentioned, the applications were written 
in Mathcad 3.1 due to sponsor requirements. Mathcad 3.1 is 
compatible with Mathcad 4.0, but Mathcad 4.0 is not compatible 
with Mathcad 3.1. 

For this thesis, a number of the antennas could not be 
adequately researched through the use of texts. In some 
cases, such as the caged dipole, very little information on 
the antenna could be found in either texts or in professional 
journals. For the caged dipole, the antenna is modeled as a 
Circular array of center-fed, equally excited, uniformly 
spaced, thin-wire dipoles to take into account the number of 
conductors that are present in the caged dipole. This model 
provides good results when compared to predicted (the caged 
dipole’s performance is Similar to that of a thick cylindrical 
dipole) and ELNEC results. 

It was found that Mathcad is not the ideal tool to 
analyze antennas due to its long processing time. An 
additional drawback of Mathcad is that when a function is in 
the beginning of an application and the same function is 
placed later in the application, Mathcad does not store the 
values of the function, but re-evaluates the function. This 
re-evaluation increases processing time significantly, 
especially when the function itself takes a very long time 


(hours vice minutes) to reach a numerical solution. 


100 


Nevertheless, the Mathcad applications developed in this 
research allow the user to analyze several antenna types 
(i.e., arrays with isotropic point sources, parabolic 
reflectors with the various feed antennas) that could not 
easily be evaluated using some antenna numerical analysis 
programs such as NEC, ELNEC, or MININEC. For the other 
antenna types researched, the antenna numerical analysis 
programs mentioned would be easier to setup, run, and provide 
parameters closer to those of a real _ antenna. The 
Gisadvantage of this approach is that it requires an ability 
to evaluate antennas using a program such as ELNEC, while the 
Mathcad applications do not require any advanced programming 


knowledge on the part of the user. 
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LINEAR ARRAY 
MATHCAD SOFTWARE ARRAY LN.MCD 


This linear array application is based on N isotropic point sources 
uniformly spaced and aligned on the z-axis. Mutual interference between 
adjacent sources are ignored when calculating antenna parameters. Since 
isotropic point sources are used as the radiating elements, polaritzation 
of the antenna can not be determined. The gain can only be idealized 
when it is assumed that the antenna input resistance is equivalent to 
the radiation resistance. Bandwith calculations are not carried out as 
a result of the narrowband characteristics of arrays, and in larogewpame 
to the dependance inter-element spacing (d) has with the wavelength (A). 
The following limits should be adhered to when analyzing Broadside, 
End-Fire, and Phased (scanning) antennas: 


0 <6, < ma (Scanning Antennas) 


@, = 1/2 (Broadside Antennas) 


05 O or n (Ordinary End-Fire Antennas) 
(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
parameters do not contain SUDSCripl ce. co, will be written AFn).) 


The linear array antenna Mathcad application will compute the following 
parameters (Items with * indicate parameters that can be CoOmpulbecmuE 
their results are idealized as the result of Using TSotropic pom 
Sources as Che wradiat tice -emeneoue 


kK = Wavelength 

Emin ~~ Minimum Distance to vene ws ar—t rele 
B = Phase Difference Between Adjacent Elements 
w = Array Factor Phase Shift 

AF, = Normalized Array Factor 

U(@) = Radiation Intensity 

Praq = Radiated Power 

Dow  EEsecelvasey 

EIRP = Effective Isotropic Radiated Power 
R, = Radiation Resistance 

Rin = Antenna Pnpueeresistanccoums 

[ = Voltage Reflection Coefficient * 

fe = TANCenNaw rte elena aa, 

G =) Gain ~ 

Ae, = Maximum Effective Aperture * 

Nem = Maximum Effective Height * 
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The following known or estimated data must be entered: 


N = Number of Isotropic Radiating Elements 

fmmomerrequency of Interest 

d = Distance Between Adjacent Elements 

fae Direction of Main Lobe 

eee antenna Feed Current 

Wee Characteristic Feed Impedence 
PRecoEelopuE data were 

N ‘= 10 (elements) d-= 0.25 

: 8 ; 

f:=3-10 CH Zz) fo := 0 

Io := 1.0 (A) Z0-= 75.0 


Define constants and calculate wavelength: 


¢ i= 2.9979-10° (meters/sec) no := 120-1 
aC PT 
Ase (meters/cycle) k= 
f A 
A = 0.999 (meters/cycle) ecd,-= || 


203 


(meters) 


(radians) 


(02) 


(dimensionless) 


Calculate linear array parameters: 


Define angularsetisceuiua 


a:=(0- 104 ae 164 


Minimum Distance to Far-Field rv:-: 


Se 


tao, (meters) 
r,:=5-N-d (meters) 
2:(N-d)? 
(a= (meters) 
Z 
He 
rmin := max(r) (meters) 
rmin = 12.509 (meters) 


Array Factor AF: 


B :=-k-d-cos( 60) 


w(8) := k-d-cos(@) + B 


snl S40) 
AFn(@) a = Ne 


snl v) 
2 
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(radians) 


(radians) 


(radians) 


(dimensionless) 


Peeration Intensity U: 


U(6) := (AFn(6))* (W/solid angle) 


Radiated Power Prad: 


7 

Prad ‘= 2-1: U(6)-sin(@) dé (W) 
“0 

Prad = 1.256 (W) 


bamectivity Do: 


Umax -= | (W/solid angle) 

Do = Am Umax (dimensionless) 
Prad 

Do = 10.007 (dimensionless) 


Effective sot ropiic Radiated Power EIRP: 


EIRP := Prad-Do (W) 


EIRP = 12.566 (W) 
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Radiation Resistance 5 





Rae — 
(Tol ) 
Rr = 2.512 


Input Resistance Rini 


Rin := Rr 


Rin = 2.512 


Voltage Re lection Gocttecerent a = 


rcs Rin - Zo 
Rin + Zo 
lT =-0.935 


Retlection EGrrcrency ean. 


a zeed {1 (Ir])"| 


a—0 125 
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(02) 


(Q2) 


(Q2) 


(Q2) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Gain G: 


G-=et:-Do (dimensionless) GdB .= 10-log(G) 


i= 1,255 (dimensionless) GdB = 0.986 


Maximum Effective Aperture Rem: 


EG 


4-n 


Aem := 





Aem = 0.1 





Maximum Effective Height Nem: 
Rr-Aem 
hem :=2- |———— 
no 


hem = 0.052 


OT 


(dB) 


(aB) 


(m) 


NORMALIZED ARRAY FACTOR RECTANGULAR PLOT 
Number Of radiating elements: 
N=10 
Direction Of Main Gobe (radians 


Jo = O 





RRR 

BERREEEDY 
mb! 

we 


0 0.628 1.257 1.885 2.513 3.142 377, 4.398 5.027 5.655 6.283 


0.75 






Larn 6)] 0.5 


0.25 
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NORMALIZED ARRAY FACTOR POLAR PLOT 





Number of Elements: =10 
Direction of Main Beam (radians): Go = 0 
, 1 ’ 1 
X(0) := |AFn(6)| -cos O+— ¥(6) := |AFn(6)|-sin[0 + - 


ec ae 


Cage Se ee 


at 
Me es. 
¢ 
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PLANAR ARRAY 
MATHCAD SOFTWARE ARRAY PL.MCD 


This planar array application is based on M and N isotropic poine 
sources uniformly spaced and aligned on the x-axis and y-axis 
respectively. Mutual interference between adjacent elements are 
neglected when calculating antenna parameters. Since isotropic pomae 
sources are used as the radiating elements, polarization of the antenna 
can not be determined. The gain can only be idealized when it is 
assumed that the antenna input resistance is equivalent to the radiation 
resistance and the total antenna efficiency is calculated for a ideal 
lossless antenna. 


(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
parameters do not contain subscripts (i.¢€., AF, will be written sane 


The planar array antenna Mathcad application will compute the following’ 
parameters (Items with * indicate parameters that can be computedmamn 
their results are idealized as the result of fusing 1sotroprc Comme 
SOUrces aS Ene radgtade mo etetemeoi: 


hk = Wavelength 

L, = Antenna Dimensiom Inyx-direct ron 

Ly = Antenna Dimension in y-direction 

Imin = Minimum Distance to the Far-Field 

B, = Progressive Phase Shift Between Adjacent Elements in x-direction 
By = Progressive Phase Shift Between Adjacent Elements in y-direction 
Wy = Array Factor Phase Shi tte ou rection 

Wy = Array Factor Phase Shi fttyin sci receion 

AF, (0,9) = Normalized Array Factor 

U(86) = Radiation Intensity 

Praqd = Radiated Power 

De = ee De reGtal yey 

EIRP = Effective Isotropic Radiated Power 

Re =.» Radgtation Resistance 

R;, = Antenna Input Resistance * 

[ = Voltage Reflection Coefficient * 

ee = “Antenna BLiverency  s 

G = Gain * 

Rem = Maximum Effective Aperture * 

Nem = Maximum Effective Height * 
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The following known or estimated data must be entered: 


M = Number of Isotropic Radiating Elements in x-direction 
N = Number of Isotropic Radiating Elements in y-direction 
f 
d 


Frequency of Interest 
= Distance Between Adjacent Elements in x-direction 


x 
dy = Distance Between Adjacent Elements in y-direction 
0, = Direction of Main Beam at @ = @o 
Seeoeerrection of Main Beam at $= 90 
ieee onctenna Feed Current 
Zo -— Characteristic Feed Impedence 
7 = Number of Increments in Degrees for Far-Field Radiation Pattern 


Enter input data here 


M :=5 (elements) dx :=0.5 (meters) 
N =5 (elements) dy =0.5 (meters) 
f=310° (Hz) Go = = (radians) 
lo = 1.0 (A) to = 7 (radians) 
Zo :=75.0 (02) j:=180 (degrees) 
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Define constants and calculate wavelengtn: 


c := 2.9979-10° 
i= g 

f 
1 = 0.999 


Calculate linear array parameters: 


(meters/sec) 


(meters/cycle) 


(meters/cycle) 


Define angular offset @ andg: 


@ =(o- 104, 10 2 aC ace 


b:=(0- 104, [0 ‘ 


J 


J 


22). Gee) 


Minimum Distance to Far-Field Cmins 


r= 1.6-4 


(meters) 


ris 5a (M-dx)? + (N- dy)? 


rmin ‘= max(r) 


2fdoan? rer i 
d 


(meters) 


(meters) 


(meters) 


Jail 


no := 120-1 
pee 20 

d 
ecd ‘= | 
rmin = 25.018 


(02) 


(m71) 


(dimensionless) 


(radians) 


(radians) 


(meters) 


Array Pacclor AFy: 


Bx :=-k-dx-sin(8@0)-cos( $0) 
By “=-k-dy-sin(60)- sin($o) 
wx(6,) :=k-dx-sin(@)-cos($) + Bx 


wy(8,) “=k-dy- sin(@)-sin(o) + By 


snl yc, sin 0.4) 


M:sin a N- sin jee) 
Z 2 


AFn(8,6) := 


Radiation imiensitey U: 


U(8, 9) = AFn(6, >): AFn(8, 9) 


Radiated Power Prad: 


| 
2 
Prad := U(6,6)-sin(6) dO do 
“O 0 
Prad = 0.206 


(radians) 


(radians) 


(radians) 


(radians) 


(dimensionless) 


(W/solid angle) 
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DiRece vine, Dae 


Umax :=1 
fees 4-n Umax 

Prad 
Do =61.09 


Effective TSor rope Radiated Power EIRP: 


EIRP = Prad-Do 


EIRP = 12.566 


Radiation Resistance Ry: 





ee stata 
( [Io] ) 
Rr =0.411 


INnpUEeReslocance Rin: 


Rin ‘= Rr 


Rin = 0.411 
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(W/solid angle) 


(dimensionless) 


(dimensionless) 


(W) 


(W) 


(Q) 


(02) 


(Q) 


(Q) 


fmteeage Relection Coeffecient [: 





i -Rin- Zo (dimensionless) 
Rin + Zo 
l =-0.989 (dimensionless) 
Peevecction Efficiency ¢e;: 
2 | 

a seed (1- (|Fl)7] (dimensionless) 
et = 0.022 (dimensionless) 
Gain G: 
G :=et-Do (dimensionless) GdB ‘= 10-log(G) (dB) 
G = 1.326 (dimensionless) Gd = | 225 (dB) 
Maximum Effective Aperture A.y: 

Ge 
Aem ‘= (m2 ) Aem =0.105 (m ) 

4-n 





Maximum Effective Height Nem: 


hem ‘= 2. aso (m) hem = 0.021 (m) 
70 


EN tae: 


TWO-DIMENSIONAL ANTENNA PATTERN: (x-z) plane 


Number of Elements along x-axis: 
Number of Elements along y-axis: 


Direction of Main Beam: 


X1(6) ‘= |AFn(0,0)|-cos ar 


cm) 
' 
t 
' 
t 
' 
' 
ee | Se 
! 
t 
] 


es nh eee 


‘ 
' 
‘ 
‘ 
' 
' 
‘ 
‘ 
ry 
‘ 
‘ 
¢ 
¢ 
' 
‘ 
t 


M =5 
N =5 
Bo = 0.524 (radians) 
do = 0.785 (radians) 


Y1(6)*=|AFn(0,0)|-sin 9-2 





xref, X1(8) 


TWO-DIMENSIONAL ANTENNA PATTERN: (y-z) plane 


Number of Elements along x-axis: M=5 
Number of Elements along y-axis: N =5 
Direction of Main Beam: Go = 0.524 (radians) 
go =0.785 (radians) 
Tl Tl Tl Tl 
200) = |AFni6,—!|-cos|6— — Y2(6) =|AFn/6,—/]-sin/@— — 
2 2 Z 2 
6 = 0° 


a as 


Te NaN a ll a allie (hls Sie 





=a | 0 ] 
xref. X2(9) 


2 us 


TWO-DIMENSIONAL ANTENNA PATTERN: $ = 1/4 





Number of Elements along x-axis. M =5 
Number of Elements along y-axis: N=5 
Direction of Main Beam: 00 = 0.524 (radians) 
do = 0.785 (radians) 
X3(0) := |AFn{@,—] | -cos|@- — Y3(0) := |AFn|6,—||-sinfo- — 
4 2 4 2 
0 = 02 


Ne | 


ew = f®@eeeewy 


a+? 


dewo2 = 


eeeeeeererereenwreeeenrereeeecercoe 





= 5 
xref, X3( 8) 
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TWO-DIMENSIONAL ANTENNA PATTERN: ¢ = bo 


Number of elements along x-axis: M=5 
Number of Elements along y-axis: N=5 
Direction of Main Beam: §o = 0.524 (radians) 
do =0.785 (radians) 
X4(6) := |AFn(6,0)| -cos Oat Y4(0) := |AFn(6,$0)|-sin - = 
pr= 302 


Se Greate a he kee ee eee eae eo ee 


ssletstiatioattal ulteineineiatios 


i] 
* 


ae 


aes 


( 
Av 


—_ 





-} 0 1 
xref, X4( 9) 
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CIRCULAR ARRAY 
MATHCAD SOFTWARE ARRAY CI.MCD 


This circular array application is based on N isotropic point sources 
uniformly spaced on the x-y plane along a circular ring of constant 
radius. Mutual interference between adjacent elements are neglected 
when calculating antenna parameters. Since isotropic point sources are 
used as the radiating elements, polarization of the antenna can not be 
determined. The gain can only be idealized when it is assumed that the 
antenna input resistance is equivalent to the radiation resistance and 
the total antenna efficiency is calculated for a ideal lossless antenna. 


(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
parameters do not contain subscripts (ite., AF, will be writteneeaaye 


The circular array antenna Mathcad application will compute thefoiic mag 
parameters (Items with * indicate parameters that can Decomputed game 
their results are idealized as the Fesult Gems ing tsorrep cero as 
SOURCES (a5 the radiating 6 lementan. 


kA = Wavelength 

Imin = Minimum Distance Co thesrar raceme 
a@, = Phase Excitation of thesmen te kemene 
AF (8,6) = Normalized Array Factor 

U(9,6) = Radiation Intensity 

Praq = Radiated Power 

Doms) DL ECEE ay ey 

EIRP = Effective Isotropic Radiated Power 
Ry = Radiation Reststanee 

Rin = Antenna Input Resistance * 

[ = Voltage Reflection Coefficient * 

f& = Antenna Efticrency.= 

G = Gain * 

Rem = Maximum Effective Aperture * 

Nem = Maximum Effective Height * 
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The following known or eStimated data must be entered: 


N Number of Isotropic Radiating Elements 
f = Frequency of Interest 
a = Radius of circle 


0, = Direction of Main Beam at 6 = @o 

Soe) Direction of Main Beam at § = $0 

ieee enicenna Feed Current 

Weeeeeeciaracteristic Feed Impedence 

feeo—menumber Of Increments in Degrees for Far-Field Radiation Pattern 


Enter input data here 


N := 10 (elements) pe (meters) 
2:1 
f:=3-10° (Hz) Qo := 0 (radians) 
Io ‘= 1.0 (A) go := 0 (radians) 
Zoe= 15.0 (Q) j:= 360 (degrees) 
Beene Constants and calculate wavelength: 
¢ = 2.9979-10° (meters/sec) no := 120-1 (Q) 
ae ; 2-0 na | 
a imerers/ cycle) k = — (m= ~) 
ry 
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Define constants and calculate wavelength: 


4 =0.999 (meters/cycle) ecd ‘= | 


Calculate linear array parameters: 


Define angularsottseuaurancn. 


g:=(o- 10°), one 22). Gee WOR) 


6 = (o- o,f 1044 24). Re iG 


Minimum Distance to Far-Field Ymin: 


- =1.6-A (meters) is =5-(2-a) 
Ps 
ae (meters) 
2 
d 
rmin ‘= max(r) (meters) rmin = 20.278 


Array Factor AF: 


fie= 12 N (elements) ie are Ray ea — 


Ca k-a-sin(00)-cos{ go - o } 
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(dimensionless) 


(radians) 


(radians) 


(meters) 


(meters) 


(radians) 


(radians) 


BeLay Factor AF: 


ie! J {k-a-sin( @)-cos (6 ~ Op) + o,) 
AF(0,6) = Se 


n 


Reatation Intensity U: 





U(@, 9) = AF(8, 4): AFCO, >) 


Radiated Power P rad: 


2% | — 
Prad ‘= U(6,)-sin(@) dé do 
0 40 
Prad =0.536 


Peeectivity Do: 


Umax := | 


_ 4-7 Umax 
Prad 


Do. 


Do = 23.427 


Effective PSORrOpaCc Radiated Power EIRP: 


EIRP := Prad-Do 


BIRP = 12.566 
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(dimensionless) 


(W/solid angle) 


(W) 


(W) 


(W/solid angle) 


(dimensionless) 


(dimensionless) 


(W) 


(W) 


Radiation Resistance Rr: 





Re ae 
( [Io] ) 
Rr— 1073 


Input Resistance Rin? 


Rin := Rr 


Rin = 1.073 


Voltage Relection Coeffecient i: 


r:- Rin- Zo 
Rin + Zo 
lf =-0:972 


Reflection  biererency (f+ 


a zed 1—irly7] 


et = 0.056 
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(Q) 


(Q) 


(22) 


(Q) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Gain G: 


Ge=et:Do (dimensionless) 


G = 1.303 (dimensionless) 


Maximum Effective Boerture Rem: 


6.15 


4-n 


Aem ‘= (m2) 








Maximum Effective Height Dem: 


hem .= 2: eae (m) 
no 





iK2 5 


GdB -= 10-log(G) 


GdB = 1.149 
Aem =0.104 
hem = 0.034 


(dB) 


(AaB) 


(m=) 


(m) 


TWO-DIMENSIONAL ANTENNA PATTERN: (x-z) plane 


Number of Elements in x-y plane: 


Direction of Main Beam: 


X1(6) = [AF(0,0)|-cos is 


¢ 
U 


0 
xref, X1(8) 
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N =10 
6o =0 (radians) 
0 =0 (radians) 


Y1(0) := |AF(0,0)| sin a 


% 


oe (Re 





TWO-DIMENSIONAL ANTENNA PATTERN: (y-z) plane 


Number of Elements in x-y plane: N =10 
Direction of Main Beam: 80 = 0 (radians) 
go =0 (radians) 
X2(8) = are -cos[6— — Y2(6) =|AF/6,—||-sin[e- — 
2 2 2 y 
8 = 0° 





6 
Qle-----= =e ee a = = Se SS y Oe 
‘7 : 
§ 
L) 


t 
4 
a 
4 


xref .X2(8) 


ie / 





TWO-DIMENSIONAL ANTENNA PATTERN: $ 





San Ue! 
Number of Elements in x-y plane: N = 10 
Direction of Main Beam: 6o = 0 (radians) 
do = 0 (radians) 
X3(0) = |AFI0,—|] -cos[o- — v3.0 :=|ar{o2 -sin{Q- — 
Zz 4 2 
@ = 0° 


eeere F-=xew# ef "=f" ffr = 


i] 
4 
i] 
i] 
6 
ay 
? 
§ 
? 
4 
? 


0 


xref, X3( 6) 
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TWO-DIMENSIONAL ANTENNA PATTERN: $ = 6, 


Number of Elements in x-y plane: N =10 
Direction of Main Beam: Go =0 (radians) 
go = 0 (radians) 
X4(6) := |AF(0,40)| -cos| @- : Y4(0) := |AF(0,40)| -sinjo- = 
pe 0- 





xref x4(8) 
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FOLDED DIPOLE 
MATHCAD SOFTWARE FOLDED.MCD 


To provide impedance matching to coaxial input cables and maintain good 
directional pattern characteristics, a folded dipole is used instead of 
a Single dipole element. The geometry of a folded dipole is such that 
the spacing between dipole elements (d) is much less than the wavelength 
(A) and the length of an element (L). 


This Mathcad application calculates the parameters of a folded dipole 
aligned vertically on the z-axis and a folded dipole positioned 
horizontally above the earth. For the calculations, the antenna is 
assumed to behave as a cylindrical dipole with equivalent radiusmiaee 


(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
parameters do not contain subscripts (1.e., Rj, will be written shames 


The folded dipole antenna Mathcad application will compute the following 
parameters {(* indicates antenna parameters that are calculated for a 
folded dipole iti beemsOace mac positioned horizontally over thewvearune 


4 = Wavelength 
— = Equivalent Radius 
Rr = Radiation Resistance 
Ra = Input Resistance of Lineareeipore 
XX, = Reactance of Linear Dipole 
Xq = Input Reactance of Linearepepole 
Zqg = input Impedance of Lineargpuagole 
Zo = Characteristic Impédance ter lo wr cell anctioslcu sain] 
Z; = Characteristic Impedance Transform of Folded Dipole 
4;n = Characteristic Impedanceme terol dedebipele 
Zg = Maximum Current of Folided wea pore 
Imin = Minimum Distance to the Far—rield 
E(6) = Electric Field * 
E,(8) = Normalized Electric Field * 
U(@) = Radiation Intensity * 
Umax = Maximum Radiation Intensity * 
Fad  — Radiaced bower 7 
Do = .Ditectivie 
EIRP = Effective Isotropic Radiated Power * 
Rin = Antenna Input Resistance * 
[ = Voltage Reflection Coefficient * 


0 


Paemerencenna Efficiency * 


Sa. Antenna Polarization Vector * 

PLE = Polarization Loss Factor * 

G = Gain * 

Aen Maximum Effective Aperture * 

eae Maximum Effective Height * 

Erp = Relative Complex Effective Dielectric Constant of Ground 
wae ertical Reflection Coefficient of Ground Plane 

Ue = Horizontal Reflection Coefficient of Ground Plane 


The following known or estimated data must be entered: 


f = Frequency of Interest 

a Radius of Conductor 

Pee RaAgiusS of Conductor 

qd = Distance Between Conductors 

L = Length of Antenna 

h = Height of Antenna Above Earth 

Peemeeeconduction/Dielectric Efficiency 

feeeerkelative Dielectic Constant of Ground 

wee eecondauctivity of Ground 

aaa far Field Distance 

1 = Number of Increments in Degrees for Far Field Radiation Pattern 
—eeeiicoming Wave Electric Field Unit Vector For Free Space Antenna 
aera Ncoming Wave Electric Field Unit Vector For Horizontal Antenna 
Ooh = Antenna Polarization Direction For Horizontal Antenna 
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ow := 10 
0 
f:=3.10° 
a := 0.0005 
b ‘=a 
d:= 12.5-a 
E=05 
h:=.5 


Enter input data here 


(dimensionless) 


(Hz) 


(meters) 


(meters) 


(meters) 


(meters) 


(meters) 


0 
owh ‘= 1/0 


rff = 1-10° 


(ca= 1 


(dimensionless) 


(radians) 


(meters) 


(radians) 


(dimensionless) 


(siemens/meter) 


(dimensionless) 


Define constants: 


oa -=|0 (dimensionless) 

0 
cis 2.9979. 10° (meters/sec) 
ee (meters/cycle) 
Lee ee, (meters/cycle) 

| -9 

for= ——-- |0 (E/m) 

36-2 


Bemepe angular offset @: 


g:20,2% 21 
i 


Meerme Cosine [Ci(x)] and sine [S1i(x)] integrals: 


n:=1]..50 


ie ] a Le 
(2-n)-(2-n)! 
n 


Sie 0.57721 + In(x) - 


no ‘= 120-1 
cet 

A 
y :=0.57721 
z-=0..2+1 


Si(x) -= 


333 


(02) 


(dimensionless) 


(radians) 


(radians) 


i= yea 1) 


ne lyon] 1) 
n 


FOLDED DIPOLE ALIGNED ON THE Z-AXIS IN FREE SPACE: 


Eguivalent radiuiswacd 


= 


a’.In(a) + b’-In(b) + 2-a-b-In(d) 


ae = exp (meters: 
(a+ b)* 
ae =0.002 (meters) 
Elect rice ire tier 0): 
_j ckerft cos{E costo) - cos EE 
E(@) =; io ae ee re (V/m) 
2-a-rtf sin( 9) 
En(6) aw, (V/m) 
2 
Minimum Distance to Far-Field Imin: 
rela (meters) aS eHle (meters) 


e206 i 
SS (meters) 
d 
rmin ‘= max(r) (meters) 
rmin = 2.5 (meters) 
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foe) 


Radiation iiecensity U| 





U(@) := fE(@)| ) (W/solid angle) 


2-no 


Radiated Power Brad 


I 
Prad := 2| U(@)-sin( 8) dé (W) 
!0 


Prad = 36.64 (W) 


Preece ivity Do: 


unas = Uf (W/solid angle) 

2 

Do mee nax (dimensionless) 
Prad 

ere 1 641 (dimensionless) 


Effective Isotropic Radiated Power EIRP: 


EIRP := Prad:-Do (W) 


EIR —60.132 (W) 
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Radiation Resistance Ry: 


Rr ene Ci(kL) + —-sin(k L)(Si(2kL)- 2-Si(k-L)) ... 


+4 cos(k-L): (1 of + Ci(2:k-L) - 2-Ci(k-L) 
2 2 


Rr = 73.281 


Input Radiation Resistance of linear qioore Rg: 


° 
° 
— 


Rd ‘= Ss 
{S) 
sin| —— 

2 

Rd = 73.281 


Reactance of linear aaipole Xi 


Xr 222 .[2.Si(k-L) 4 cos(k-L)-(2-Si(kK-L) — Si(2-K:L)) .. 
4-n 





Ar = 43.142 


input reactance of linear dipole Xj: 


° 
. 
— 


Ge 
( | call 
smnmj-— 
5 
Xd = 43.142 


1S '6 


(02) 


(02) 


(02) 


(02) 


(02) 


(Q) 


(02) 


(02) 


Input impedance Se Neatealpole Zqi 


. 
* 
—— 


Zd =Rd4j -Xd (Q2) 


Zd = 73.281 +43.1421 (02) 


Characteristic impedance of two-wire transmission line Lie 





| a | (2) 
: 2-.Ja-b 
Zo = 302.312 (02) 


Characteristic tmpecdicesCranchorm tommrolded dipole Jet 


Zt := 43 Zortn( A (02) 
° 
Zt =-2.747 10°3 (Q) 


Characteristic impedance of folded Gipole Zins 


4-Zt-Zd 
= ——___ (22) 
2-Zd + Zt 
Zin = 293.307 +172.4651 ((2) 
Input Resistance Rjy: 
Rin = Re(Zin) (02) Roi 73.507 (Q) 


5 / 


Voltage Relectlonmeoct tecteneml 


r _ Rin- Zo 
Rin + Zo 


(dimensionless) 


Reflection Efftictene see 


et = 0d-11- (Ir1)°] (dimensionless) 


Gain G: 
G:=et-Do (dimensionless) 
G = 1.64] (dimensionless) 


Polarization .0ss factor (Pilr- 


wee 1 
PLE = (low cal} (dimensionless) 


Maximum Effective Aperture Rem: 


Cie 


4-n 


Aem ‘= 





-PLF (m2 ) 





Maximum Effective Height Nem: 


Toney) Se 
nO 





P==O:0i5 


et = | 


GdB := 10-log(G) 


GdB =2.15 
lee 

Aem = 0.13 
hem = 0.318 


ree 


(dimensionless) 


(dimensionless) 


(dB) 


(daB) 


(dimensionless) 


(m@ ) 


(m) 


Far Field radiation pattern of vertical folded dipole in free space. 


X(@) = JEn(0)|-cos 0+— Y(6) := [En(@)] -sin po 


yrrort ccc coe eee eee 


wer = 


Frew 


¢ 


i] 
i] 
i] 
i] 
i] 
Cr 


eewerere =; 





ie 2 


HORIZONTALLY POSITIONED FOLDED DIPOLE OVER THE EARTH: 


Define angular offsetseymeand 102; 


(radians) 02. 2-447 


1 
Lar 
2.1 2 21 


i 
2 


Relative complex jeffectiverdielect rrc comarane Erp? 





ap = er J 
2-1-f:0 


¢ component of E-field Eg: 





sj) eit cos sinc sino) - cof 
E4(@1) :=j -no— ee ee ee ae 


2-n-rtft . ae 2 
] = sin(@1 )°-sin(0) 


kl , k-L 
cos] ——-sin[{@2_}-sin(O)} - cos|— 
ela 2 ( ) | | 2 | 


mice i = sin(02,)  sin(0)? 





Eo] =) “No 


§@ component of E-field EQ: 


eae cos} —-sin( 61! inf) = cos 
a of a ae 


E0(@1) :=j -jo<——- 
1— sin(0l sin 
z 





tert co AE sino) sof] - cos 
_ eJ 2 2 2 
E01] :=j -no eee 


i 2-n-rff 
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(radians) 


(dimensionless) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 





Wemmebeal plane wave reflection coefficient I,: 


rv(él)*= erp-cos( 81) - Alba sin(@1 y 














(aamens Hon Less, 
erp-cos( 81) + aferp - sin(@1) 
erp-cos( 02 } - = sin(o2,)” 
—————— (dimensionless) 
erp-cos{ 62 + Jjerp- sin(02,}° 
Pommzontcal plane wave reflection coefficient I): 
os(01) d n(O1)° 
pe) — wep sintlt) (dimensionless) 
cos( 61) + ee sin( 6} )? 
cos(62 } - ap sin(02,) 
——————————————————————————— (dimensionless) 
cos(02 } 1 ap - sin(o2_)° 
Peete la (6 Component) with ground plane reflection Edgp 
Edgp(01) -= 401): (ed KP). rag] yer) -Hh-eos( 80) (V/m) 
i -k-h-cos{02 Bye ich: 
ct, +n, sm, =) Oe 
Egn(o1) = JEEP (v/m) 


max(Ed¢gp 1) 
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E-field (9 component) with Ground plane (ret lec unas EQgp: 


E0gp(61) :=E0(01)-(e k-h-cos( 61) _ rv(0l)-€% kehreos(01)) 





|Edgp(61)| 
max(E@gp1 ) 


E6n(01 ) := 


Radiatien Intensity Unig. 


Uh(81 ) =| tetepc01 yy? + (AB0gp¢ 0191)" 
‘T!0 


Radiated Power PRrad: 


fu 
2 
Phrad = 2-1. Uh(61 )-sin(61 ) d61 
0 
Phrad = 70.415 


Pre criviitry Diy: 


a 
EOgp | 





(liven 


) + 











Ugp maxX... 
_ Zz 
2-no 


| 


Uhmax = max(Ugp_max) 
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(V/m) 


(V/m) 


(V/m) 


(W/solid angle) 


(W) 


(W) 


(W/solid angle) 


(W/solid angle) 


Perectivity Bee 


Dho 222 Uhmax 
Phrad 
Dho = 3.938 


Effective Isotropic Radiated Power EIRPh: 


EIRPh := Phrad: Dho 


EIRPh = 277.3 


Input Resistance RAint 


Rhin ‘22 nad 
(st) 

sin|— 

2 

Rhin = 140.83 


eeeage Relection Coeffecient [H: 


rH cs eeu ZO 
Rhin+ Zo 
TH =-0.364 


Bepbection Efficiency ehy: 


eht 200d 1 — ( [rH] | (dimensionless) eht = 0.867 
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(dimensionless) 


(dimensionless) 


(W) 


(W) 


(Q) 


(Q) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Gain Gh: 


Gh := eht-Dho (dimensionless) 


Gh =3.415 (dimensionless) 


Polarization loss factor PLFv: 





Ehmax ‘= |E@n(@ph)| (V/m) 
0 
0 . 
oah := (dimensionless) 
Ehmax 
Lobes 


PLFh ‘= (Peal 


PLFh = 7.736 10 ~ 


Maximum Effective Aperture AVem: 


Gh? 


4-n 


Ahem ‘= -PLFh (m= ) 





Maximum Effective Height Vem: 


pesee on Rhin- Ahem eat 
10 
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GhdB ‘= 10:log(Gh) (dB) 


GhdB = 5.334 (dB) 


Lobes ‘= toor[2* + } 
d 


(dimensionless) 


(dimensionless) 


Ahem = 2.099" 10 (m2) 


hhem = 0.018 (m) 


Pemlanew Radiation Pattern (y-z plane) of Horizontal Folded Dipole Over 
The Earth: 


Y1(61) = [Egn(61)|-sin{61 ae X1(01) = |E9n(61 )| -cos| 61 2s 


Fae See, ale 





xref, X1( 01) 
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E-Plane Radiation Pattern (x-z plane) of Horizontal Folded Dipole Over 
The Earth: 


Y2(01) ‘= |Een(61)|-sin{ol +s X2(61) “= |E@n(01)|-cos{ 01 + . 


Rete tate t aster tates 


Y2(91) 


oe | 





xref_,X2( 01) 
m 
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CAGED DIPOLE ANTENNA 
MATHCAD SOFTWARE CAGED DI.MCD 


ieemeagea dipole application is modeled after a circular array with the 
Poeederaig elementS as dipoles. The dipoles are modeled after conductors 
with infinitely small radii and antenna length greater than A/4. Mutual 
interference between adjacent elements are neglected when calculating 
a@mtcenna parameters. 


(Note: Mathcad equations can not use symbolic subscripts. Therefore, 
Bete cers GO not contain subscripts (1.e., AF, will be written AFn).) 


mae caged Gipole antenna Mathcad application will compute the following 
parameters: 


VERTICAL CAGED DIPOLE IN FREE SPACE: 


A = Wavelength 

eeaeeoinimum Distance to the Far-Field 
AF (66) = Normalized Array Factor 

E(06) = E-Field of Finite Length Dipole 
awe = Total E-Field 

Eee) = ‘Total Normalized E-Field 
U(86) = Radiation Intensity 

Peace Kadiated Power 

Pees Oirectivity 

BIRP = Effective Isotropic Radiated Power 
iaaeaeenadiation Resistance 
foeaeee)|«6CLMput Resistance 

Tr = Voltage Reflection Coefficient 
feaoeeontcenna Effeciency 

G = Gain 

PLF = Polarization Loss Factor 

pene) Maximum Effective Aperture 
iene) Maximum Effective Height 

flow = Low Frequency 

mega —= High Frequency 

BW = Bandwidth 


VERTICAL CAGED DIPOLE OVER EARTH: 


Erp = helm secomplex ELrective Dielectric Constant 


AFv (81,6) = Normalized Array Factor 
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Ev (61,6) = E-Field of Finite @ienigthepueel. 


r,(01) = Vertical Plane Wave@Retleemremm coer amen 
Ev ei0l ge = Totalaeeraetc 

EV, , (Ole) = eta l® Nowmerli Zecies— Fi elk 

Uv (61,6) = Radiation Intensity 

PVyrsaq = RdGiatedsrever: 

DVG° = PPLE eccer, fey 

EIRPv = Effective wiseotropic Radiorederc ie: 
Rv, “= “Radiata ongnesmeranec 

RVj, =) JNPUL Rest= tance 

[Vv = Voltage Reflection Coefficient 

eV, Samepmtennay Ski tecleney 

Gv = Gain 

PLEV =" Polarization Goss Factor 

PVen = Maximum Bi fective Apereure 

hve, = )Maximumebttective Heugie 


HORIZONTAL CAGED DIPOLE OVER EARTH: 


AFh (81,61) = Normalized Array Factor 

fy, (61) = Horizontal Plane Wave Reflection Coefficient 
E60n (861) = 6 Component of E-Field 

E@h; (891,91) = @ Component of Total H=Firetd 

E6h;y, (81,91) = 6 Component of Total Normalized E-Field 
Eém (81) = 16 Component of E-Field 

Egh; (91,61) = $¢ Component of Total E-rierd 

E¢nt, (61,91) = @ Component of Total Normalizedmrai-ad 
Uh (61,61) = Radiation Intensity 

Parsaq = Radiated Power 

Diane. Daweel ivwivey 

EIRPA = Effective Isotropic Radtatedseeve: 

Ri; - = Radiation Resistance 

Riga = lbpUG Resicrance 

[H = Voltage Reflection Coefficient 

fhe = “Alycenned EPrtretency 

Gh = Gain 

PLEA .= Polarization Loss tactem 

PNem = Maximum Effective Aperture 

hhem = Maximum Effective Height 
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The following Known or estimated data must be entered: 


N = Number of Dipoles 

a = Radius of Antenna 

L = Length of Antenna 

h = Antenna Height Above Ground 
Peep cequiency Of Interest 

1 eee Sntenna Feed Current 

Z 


° 
pea lnput Impedance 
fence lative Dielectric Constant of Ground 
foe—meeconauctivity of Ground 
eeaerncoming Wave E-Field Unit Vector for Antenna (Free Space) 
Seeeeeicenna Unit Polarization Vector (Free Space) 
coemeaeelncoming Wave E-Field Unit Vector for Vertical Antenna 
OVD = Direction of Incoming Unit Vector (Vertical Antenna) 
Tien incoming Wave E-Field Unit Vector for Horizontal Antenna 
Why =— Direction of Incoming Unit Vector (Horizontal Antenna) 
j = Number of Increments in Degrees for Far Field Radiation Pattern 
ENCer inpul Gatawhere: 
l l 
ow ‘= | 0 (dimensionless) oa ‘= |0 (dimensionless) 
0 0 
l 
. 80: 1 
ovw =]|0 (dimensionless) 8vp =—— (dimensionless) 
180 
ohw :=|0 (dimensionless) Ohp :=0 (dimensionless) 
l 
N =8 (elements) Io = 1.0 (A) 
ais (meters) Zo :=75.0 (Q) 


oO |— 
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Le=05 (meters) y:=90 (increments) 


He 025 (meters) O10 (siemens/m) 
f:=3-10° (Hz) = 15 (dimensionless) 
Define constants and calculate wavelength: 
¢ := 2.9979. 10° (meters/sec) no := 120-2 (Q) 
Paes: © 2 al 
Ace — (meters/cycle) ki=—— (ma) 
f d 
eee (A) rff = 10° (meters) 
N 
ts ee (farads/m) ecd <= | (dimensionless) 
36-1 
4 = 0.999 (meters/cycle) z-=0..25 (increments) 
Calculate caged dipole parameters parameters: 
Detine angular vobrsaenu 
2: . 
6°2 0,2" ..2-1 (radians) Reais Go: (radians) 
J J 


B50 


Minimum Distance to Far-Field Imin: 


r= 1.02 (meters) 


eae + (2-2)° (meters) 


elie) (2-a)° | 
r, —— 


(meters) 
rmin ‘= max(r) (meters) rmin = 2.795 


VERTICAL CAGED DIPOLE IN FREE SPACE: 


Array Factor AF(Q): 


m= 1..N (elements) Ofer 


| -k-a-sin( 98)-costo— © 
AF(8,6) =) line’ Ks O) ont) 


n 


Meera) Of finite dipole E(0): 


kL k-L 
cos]|-—-cos(8)} - cos] — 
an | 2 | | 7 


E(@) =} jo , 
2-nrff sin(@) 
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(meters) 


(radians) 


(dimensionless) 


OU tal 


Total E-Field Ey, (0) 4 Normaluzedse neue Erp (8) : 


Et(0,$) := E(@): AF(6, 4) (V/m) 


Etn(6, 6) .-_ _Et(8,6) (V/m) 
ei} te 
2 Z 
R 


adiation Inmtensa MUO: 


U(6, 6) 2 ({Et(0, 61 9° 


DS, 


(W/solid angle) 


Radiated Power Bac 


(W) 
U6, >)-sin(6) d6 do 





Prad =0.441 (W) 


DIReculvity Do: 


umax:=U[5. (W/solid angle) ae _ 4-7 Umax 
2 


(dimensionless) 
Prad 
Do = 1.543 (dimensionless) 
Effective isotropic Radiated) Power seine. 
EIRP := Prad-Do (W) EIRP = 0.681 (W) 


2 


Radiation Resistance Ry 


_ 2:Prad 
( In] ° 


Rr: (Q2) 





Input Resistance Rin: 


2) (Q) 


Memeage Relection Coeffecient [: 


rcs Rin - Zo 
Rin + Zo 


(dimensionless) 


bemeecti0n Efficiency ef: 


et ‘eecd-| 1 - aba 


(dimensionless) 
Gain G: 
G :=et-Do (dimensionless) 
G =1.513 (dimensionless) 


Polarization Loss Factor PLF: 


PLF “= ew-ca Ir 


(dimensionless) 


53 


Rr = 56.487 

Rin = 56.487 
r=-0.141 

et = 0.98 

GdB ‘= 10-log(G) 
GdB = 1.798 
PLF = 1 


(02) 


(dimensionless) 


(dimensionless) 


(dB) 


(dB) 


(dimensionless) 


Maximum Effective Aperture Agm: 


Gum 
Aem := _PLF (m2 ) Aem =0.12 (m2 ) 
4-n 








Maximum Effective Height ent 


hem ‘= 2: Rr-Aem (m) hem = 0.268 (m) 
; 10 


Define Frequency Relationship F(x): 





F(x) =- 5.696 10° >-x4 0.31481 (dimensionless) 


Freguency Range fjoy and fhigh: 


flow :=f- ies -f (Hz) flow = 2.056 10° (Hz) 
sete. IF | 8 
a (ro tl open Fiz) thigh = 3.944 10 (Hz) 


Bandwidth BW: 


BW =fhigh- flow (Hz) BW = 1.888*10° (Hz) 
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PeetetuD PATTERN: This far-field radiation pattern is for the vertical 
Segeercipole in free space: 





Number of Dipole Conductors: N =8 
Radius of Antenna: a=0.125 (meters) 
n | n 
X(6) -= |Etn(6,4)| -cosf@- — Y(6) -= JEtn(6,6)| -sinfo- ; 
2 
b=, 0° 


T 


i 
4 
L) 


ace <i = = 


hen 
t] 
4 
é 
t 
é 
é 
t] 
i 


4 
‘ 





xref, X( 8) 


6 = 180° 


5 


VERTICAL CAGED DIPOLE OVER EARTH: 


Define angular Of Esetwul sancuce 


6] 2 a (radians) 2 = (radians) 
22 eee = 
82 =a ees (radians) 
zZ 2+] 


Relative complex ertective dieleceriecreconseane Erp: 








sp =er— J - (dimensionless) 
2:1: f-60 
Array Factor AFv(@): 
} -k-a-sin( 81 )-coslo— © 
AFv(61,9) =S ine’ ae on(p d (dimensionless) 
n 
j -k-a-sin{62,) -cosl$2,— 
AF vgp,, =} line a ’ coslt2e ) (dimensionless) 
n 
E=pield tot irpite apo lemme 
= COs ECL) — cos eae 
Ine? Kt Z 2 
Ev(@1) =) -7o————— |_—————————qx— (V/m) 
2-nrff sin(61 ) 
Bis cos ee Q2 }| - cos Ls 
eeiree ae 2 2 
Evep,*=) -\o—_—$— = | A —__——— (V/m) 
2 arff sin(62) 
Vertical plane wave reflection coetficlenul (uum 
| 2 
[veel c= ae ENGR (dimensionless) 


erp cos( 1 ) + = - sin(6] - 


To 


Evt(61,6) =Ev(6l )-AFv(61 , $)- (J “k-h-cos( Ol) rv( 61 )-e ama 


Evigp, = 








Ever fone {e Kenvcos(O2) 


[Evt(01 , 6)| 
max( Evtgp) 


Evtn(@1,6) = 


Paemeatron Intensity Uv (81): 


Uv(01,6) := Flea 1)’ 


“0 


vst = (ive 


Radiated Power PV rad: 


Dan | = 

Pvrad ‘= Uv(61,)-sin( 61) d6l do 
0 70 

Pvrad = 0.483 


bameectivity Be. 


Uvmax := Uv( 6vp,0) 


| es, tween) 
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(dimensionless) 


(V/m) 


(V/m) 


(V/m) 


(W/solid angle) 


(W/solid angle) 


(W/solid angle) 


DEeerl Vasey DV5: 


_ 4-1 Uvmax 
pve? = ——— 
Pvrad 
Dvo = 4.415 


Effective Isotropic Radvated Pover ss hia. 


EJRPv ‘= Pvrad: Dvo 


(W) 


Radiation Resistance RVs 


__ 2:Pvrad 
({In[ )° 


Rvr 





(02) 


Pape Resistance ae 


Rvin “= 


iesth = 
sim j|}—- 
2 


(Q) 


Voltage Relection Coeffecient ves 


ry v= Rvin= Zo 
Rvin + Zo 


(dimensionless) 


Reflection Eiticilencvige. 


vt :=ecd| 1 - ([rv1)?| 


(dimensionless) 


Gane: 


EIRPv =2.133 


Rvr =61.857 


Rvin = 61.857 


TV =—0 026 


evt = 0.991 


(dimensionless) 


(dimensionless) 


(W) 


(Q2) 


(dimensionless) 


(dimensionless) 


Gain Gv: 


Gv =evt-Dvo (dimensionless) 


Gv =4.374 (dimensionless) 


Polarization Loss Factor PLFv: 


Evmax ‘= Evin(@vp,0) (V/m) 


Evmax 





lobesv 
Ova -= 
O 


0 


(dimensionless) 


2 
BEFY = New oval) (dimensionless) 


Maximum Effective mBpertute Rem: 


Gv: i 


4-1 


Avem ‘= -PLFv (m4) 





Maximum Effective Height hAVen: 


Ate = 9. Rvr:Avem on) 
no 


GvdB =10-log(Gv) (dB) 


GvdB = 6.409 (aB) 


lobesv ‘= floor 2% + } 
A 


PLFEv = 0.986 (dimensionless) 
Avem = 0.343 (me ) 
hvem = 0.474 (m) 
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E-FIELD PATTERN: This far-field radiation pattern is for the vertical 


caged dipole over earth: 
Number of Dipole Conductors: 
Radius of Antenna: 


Height of Antenna Above Earth: 


X1(61) := [Evtn(@1 ,6)|-cos/o1+— 
: 


yref 


Y1(81) 


N =8 
a=0.125 
h=0.25 


(elements) 


(meters) 


(meters) 


Y1(01) := |Evin(61,$)|-sin{ 6] ee 


1 

‘ 

(] 

1 

t 

1 

t 

i] 
rc 
t 

i] 

i] 

i] 

t 

i] 

‘ 
ld 


xref XI 61) 
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HORIZONTAL CAGED DIPOLE OVER EARTH: 


Demgeme angular offset 91 and $2: 
; ui 
gl -=O,—..2-1 (radians) 


J 


Pemayeractor AFh (01,91): 


n‘=1..N (elements) 0 ‘= 2-1—— 


yj -k-a [sin -sin( 81 )-sin(© s(91)-cos(O 
AFh(OL,91) =) Ine (in 91-sin( 01 -sin(GQ) +-c0s( 81)-cos (04) 
Nn 


AFhL, ‘= Stine -k-a (sin(92) -sin(62) -sin(9,) + cos (62) -cos(0,)) 





* 
' 


kL. kL 
a cos} —-sin(@1 )} - cos AE 
-) -k-rff ( 2 | 2 


In-e 


2-n-rff . @. 
1 — sin(@1 ) 


kL. k-L 
ae cos|-——-sin{ 62 }}| — cos} — 
In-e? ae 2 ( ) | | 


toes Te sin{ 02} : 


E0h(@1) =j -no 


E@h1, *=j -no 


iereel E-Field EOh; (01,61) : 


Eéht(@1 , 61) = E@h(61 )- AFh(@1 , 61 fe) ‘Kh-cos(O1) rv el )-e) K-h-cos( 1)) 





E@htl | ‘= [em arm, (c A-h-cos (62) : Ivi-e” sia | 
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@ component of E-Field of finite dipole E@h(@1): 


(radians) 


(dimensionless) 


(dimensionless) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


Normalized B-rield seth Se xieiol 


Eohtn(61 . 61) := Leone ob ie 


max(E@ht1 ) 


Horizontal plane wave reflection Goer ric Nemem meee 


rho) := cos(@1)- Jap - sin( 61 ;. 
cos(@1) + erp — sin(01)° 





(dimensionless) 


rn S2(02) = for~ sly 
cos(02) + Joy sin(o2)’ 


¢ component of E-Field of finite dipole Ed¢h (01): 


(dimensionless) 


pee) Baa .. 
Eoh(@1) =4 whe (i - cos KE) (V/m) 
2-nrff 2 
ed KE - 
Eohl :3j wt = coo SE} (V/m) 
: 2-n rff 2 


Total E-Field E@gh; (91,41) & Normalized E-Field E@¢h;, (81,61): 


Eoht( 61,61) := E¢h(@! )-AFh(@1 , 61): (.i Ebicost Ol) aincalpree k-h-cost 61) (V/m) 


| (V/m) 





E¢htl | = | Eohl -AFh1 - (e k-h-cos (2) n hl -e” aa 4 


Eohtn( 61 , $1 ) _ [E¢n(o1 sD] (V/m) 
max( Eght! ) 


L6Z 


Eoematron Intensity Uh(01,$1): 


8, 
Uh(01 41) = [ |Eontc01 , 411 7+ ({Eehtco1 41) 197] 
‘yO 


uo = | (ets (LED 


Radiated Power Eyes 


er | — 

Phrad .= Uh(O1 , 1 )-sin(01 ) dOl dol 
0 0 

Phrad = 1.059 


Dieectilvity Dho: 


Uhmax -= Uh(6hp,0) (W/solid angle) 


Dho = 5.051 


Effective SOE EODIC Radiated Power EIRPh: 


EIRPh := Pvrad: Dvo (W) 


Radiation Resistance Ree 


_ 2-Phrad 


= ; (Q) 
([In| ) 
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Dh _ 4-1 Uhmax 

Phrad 
EIRPh = 2.133 
Rhr = 135.58 


(W/solid angle) 


(W/solid angle) 


(W) 


(dimensionless) 


(dimensionless) 


(W) 


(Q) 


InPpue Resistance RAin: 


Rhin = — eer (Q) Phi (eos (Q) 


Voltage Rellectiome Coctieerenwa! i: 


i Rhin- Zo 
Rhin+ Zo 


TH: 


(dimensionless) TH =0.288 (dimensionless) 


Reflect ron rete keney see: 


eht = 0d 1 — ( \TH| | (dimensionless) eht = 0.917 (dimensionless) 
Gain Gh: 

Gh = cht-Dho (dimensionless) GhdB =10-log(Gh) (dB) 

Gh = 4.633 (dimensionless) GhdB = 6.659 (aB) 


Polarization Loss Factor PLFEh: 





Ehmax := Eghtn(6hp,0) (V/m) lobesh ‘= Noor(22 
d 
0 
0 
cha:= (dimensionless) 
hmax 
lobesh 
ale (dimensionless) 
PLFh = ere) ee) = (dimensionless) 
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Maximum Effective Aperture ANem: 


Ghia’ 


4-n 


Ahem := -PLFh (m@) 





Maximum Effective Helgi hhem: 


hhem := 2: - A Mate (m) 
nie) 


OS 


Ahem = 35.536 


hhem = 7.15 


(m ) 


E-PLANE RADIATION PATTERN: This far-field radVativoneparcerne ene 
(x-2Z) plane is for the horizontal caged dipole over vcarrtn. 


Number Gf -Dipole Conaucrors: 


Radius of Antenna: 


Height of Antenna Above Earth: 


X3(01) ‘= |E@htn Ol -cos|@1 4 — 


yref 


Y3(91) 


Z 


N =8 (elements) 
a=0.125 (meters) 
h=0.25 (meters) 


Y3(01) := |Edhtn al -sin 1+ — 





xref, X3( 61) 


H-PLANE RADIATION PATTERN: This far-field radiation pattern in the 
Meemplane is for the horizontal caged dipole over earth: 


Number of Dipole Conductors: N =8 (elements) 
Radius of Antenna: a=Q.125 (meters) 
Hemme Or Antenna Above Earth: h =0.25 (Meters) 
X2(61) := |Eghtn(61 ,0)| -cos{@1 + = Y2(01) := |Eghtn(@1 ,0)|-sinfel + — 
7 7 


-_ ys 


4 
i) 
i) 
i] 
i} 
i] 
t] 
i] 
i] 


Teer See ee 





xref. X2( 61) 
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THE HELICAL ANTENNA (REFLECTOR OPTION) 
MATHCAD SOFTWARE-DISH HEL.MCD 

When built to the proper specifications, the helical antenna possesses 
many qualities which make it suitable for a wide variety of 
communications applications. If the following conditions are satisfied 
the helix will exhibit a highly directional axial main lobe, low side 
lobe level, negligible mutual interference with adjacent antennas, low 
voltage standing wave ratio (VSWR), and resistive input impedance over a 
wide frequency band: 


8 < Gye eis 
Kine> 
Too ld 


(Note: A in a subscript indicates the dimension is in wavelengths. 
Mathcad equations can not use symbolic subscripts. Therefore, the symbol 
A will immediately follow the parameter in equations (i.e., CA) in lieu 

of subscripts.) 


The helical antenna Mathcad app licatien will compute the following 
Parameters . ives with * indicate Patameters@enat Tare calculated for 
both axial and Peripheral feed Geometries): 


= Circumference of Helix 


= Wavelength 
a = Pitch Angle 
Do = Directivity 
p = Relative Phase Velocity 
E94 = Electric Field Components 
Y = Array Factor Phase Shift 
U = Radiation Intensity 
Prad = Radiated Power 
R = Antenna Input Resistance* 
[ = Voltage Reflection Coefficient* 
€. “= Reflect von] Bette cicy. 
he, = Maximum Effective Heignur 
GP Gane 
EIRP = “Effective D@sotrepic ad tacecueone um. 
Fey = Maximum Ef fective Aperemnes 
AR = Axial Ratio 


PLF = Polarization Loss Factor 
BW = Bandwidth 

fhigh = Upper Frequency Limit 
Liow — LOwWwer Frequency iimit 


Acceptable Conductor Diameter 
Ex, y,Z = Electric Field Cartesian Coordinates 


aG'8 


O57 bp Zoi temaolaraizatiton vector Coordinate Angles 
feeeeeencenna Unit Polarization Vector 


Imin Mininum Distance to the Far-Field 


The following data must be input based on known or estimated data: 


D = Diameter of Helix (Center to Center) 
S = Spacing Between Turns (Center to Center) 
Beaerength Along Conductor of One Turn 
n = Number of Turns 
d = Diameter of Helical Conductor 
f = Frequency of Interest 
m = Desired Mode 
1 = Number of Increments in Degrees for Far Field Radiation Pattern 
feeeeeoncenna Feed Current 
Weeeeeechiaracteristic Feed Impedance 
aeeerncoming Wave Electric Field Unit Vector 
Engel IJnpuemaatca here: 
| 
AD 
a. (dimensionless) 

ow =] J 

iF 

0 

me). ! 

D ‘= — (meters) L := 0.097 (meters) 

nt 

: ; 10 

n:=10 (turns) f := 3-10 (Hz) 
m ‘= | (dimensionless) Za = 150 (Q) 
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d =0.005 
x= ] 
y=l 


(meters) 


(meters) 


(meters) 


(meters) 


1 = 360 (degrees) 
Io =] (A) 
z := 1000 (meters) 


Calculate helical Geometric Parameters and define constants: 


¢ ‘= 2.9979-10° 
Cs) 
C=0.] 
C) £ 

r 
Ci = 10.007 
ee =e 

d 
Li = 9.70679 

d eee 

Tt 

ad = 12.95276 


(meters/sec) 


(meters) 


(meters) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(degrees) 


(degrees) 


aie 120-2 
yee he 

f 
= 9.993°10_ 


ma 
a := atan|— 
é 


a = 0.22607 


SA = 2.30161 


7 O 


(Q2) 


(meters/cycle) 


: (meters/cycle) 


(radians) 


(radians) 


(dimensionless) 


(dimensionless) 


(m7 1) 


Calculate helical antenna palbameters: 


Weermeceanguiar offset 6 Minimum Distance to the Far-Field 
from helical axis: ea 
r, eet (meters) 
eo". 2.x (radians) 
1 r i= Sn (meters) 
Tl r _ 2(n8)’ (meters) 
p20,[%].24 (radians) p aie \ 
1 
rmin ‘= max(r) (meters) 
rmin = 10.5874] (meters) 
Relative Phase VeTOeKL y. p= 
. La 
_* (dimensionless) 
See im+ |— 
n 
p =2.89616 (dimensionless) 
meeayeractoOr Phase Shift w. 
Li . 
w(@) t= 2-1 1SA-cos(6) — — (radians) 
Pp 


eat 


Electric Field Components Bg, Eg: 





sol si 26) 
2: Z 
E(6) := -cos( 8) (V/m) 
sin ¥(8) 
2 
E@(@) := E(6) (V/m) 
Fo(0) :=j -E(@) (V/m) 
Radiation “Wimeensitey Ue 
U(6) ae |E(6)| ‘c (W/solid angle) | 
No 
Radiated Power P rad: 
2n [2 
Prad ‘= U(6)-sin(@) d@ do (W) 
a) Sa 
Prad = 1.67351°10- (W) 
Di recta rey. Do: 
Do-= 12-Ch-n Sh (dimensionless) Do2 ee) (dimensionless) 
Prad 
Do =2.7658 10° (dimensionless) Do2 = 19.91822 (dimensionless) 


daz 


Axlal Ratio AR: 


ar:=|1x(sina)-“} 
Pp 


AR = 1.17586 


Effective Isotropic Radiated Power EIRP: 


EIRP := Prad:-Do (W) 
EIRP = 46.286 (W) 


Polarization Loss Factor PLF: 


| 2 2 
8p =a 47 (radians) 


Zz 


6p = 1.41421 10° 


(radians) 


Ex =E0(8p)-cos( @p)-cos( dp) — Ed( 6p )-sin( dp) 
Ey =E6(@p)-cos(@p)-sin(dp) + Ed( 8p )-cos( dp) 


Ez =E@(@p)-sin(8p)-- 1 


Ex 
oa -__. Se ee Ey 


Aclexl)?+ ClEyI)24 (IEzl 2 Ez 


05- 0.5) 
pe 0.5 + 0.5) 


-9.99999" 10 * 


EIJRP2 := Prad:-Do2 


Ne 2) =O 555 


op = anf 


x 


bp =0.7854 
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(dimensionless) 


(dimensionless) 


(W) 


(W) 


(radians) 


(radians) 


(V/m) 


(V/m) 


(V/m) 


(dimensionless) 


(dimensionless) 


Polarization Loss Factor PLF: 


—} 2 
PLE = (low-cal } (dimensionless) 


Radiation Resistance Rye 


Prad 
(|Io] )? 


Rr :=2-: (02) 





Dual Parameters 


Axial Feed 


Input Resistance R: 


Ra = 140-fC2 (Q) 





Ra = 4.42874 107 (Q) 


Votcage Re tlect ton Cocttectenum: 





Ra— ae 
[as= (dimensionless) 
Ra+ Z 
0 
Ta = 0.49399 (dimensionless) 


Retleceiron Priuevency ¢y 


sra:= 1 - ({Taly? (dimensionless) 


era = 0.75597 (dimensionless) 
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ele == | 


(dimensionless) 


Rr = 3.34702 10 ° (Q) 


Peripheral Feed 








Rp ._ 10 
Jen 

Rp = 47.41756 
Rp- Z 

Ip = 2 
Rp+ Z, 

Pp =-0.51962 


ap :=1-((pl)° 


erp = 0.72999 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Gain G: 


Ga := era Do 


Gadb := 10: log(era: Do) 


Ga = 2.09088- 10° 


Gadb = 43.20328 


Ga2 := era Do2 


Gadb2 := 10-log( era: Do2 ) 


Graz = 15.05/65 


Gadb2 = 11.77757 


(dimensionless) 


(dB) 


(dimensionless) 


(AaB) 


(dimensionless) 


(dB) 


(dimensionless) 


(dB) 


Maximum Effective Aperture Rem: 


2 
conn = LST 
4-n 
Aema =0.16615 
era-1°-Do2 
Aema2 := —————— :PLF 
4-7 


Aema2 = 1.19657 10° 
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Gp ‘= erp: Do 


Gpdb := 10:log( ap: Do) 


Gp = 2.01902"10° 


Gpdb = 43.0514 


Gp2 := erp Do2 


(dimensionless) 


(dB) 


(dimensionless) 


(dB) 


(dimensionless) 


Gpdb2 := 10-log(erp:-Do2) (dB) 


Gp2 = 1454017 


Gpdb2 = 11.62569 


2 
Aemp ‘= erp:-4 Do -PLF 
4-7 
Aemp = 0.16044 
4?.Do2 
Aemp2 sb 4 =°* pyr 
4-n 


Aemp2 = 1.15545 10 — 


(dimensionless) 


(dB) 


Maximum Effective Height Nem: 

















hema := 2: pgs (m) hemp := 2- {Rr- ally 

"lo No 
hema = 2.42912*10 * (m) hemp = 2.38701*10 > 
hema2 ‘= 2- |Rr- ail (m) hemp2 :=2- {Rr- alll 

No No 
hema2 =6.51873 10 > ” (m) hemp2 =6.40574 10 ° 
Bandwidth: 
fhigh := mee (Hz) flow cee 
C 

fhigh = 3.44759 10° (Hz) flow = 2.39832* 10° 
BW ‘= fhigh - flow (Hz) 
BW = 1.04927 10° (Hz) 


Acceptable Conductor  Prameter: 


dmin ‘= .005-A (m) dmax := .05-) 


dmin = 4.9965 10> (m) dmax = 4.9965 10 ~ 
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(m) 


(m) 


(m) 


(m) 


(H2) 


(Hz) 


(m) 


(m) 


HELICAL ANTENNA FAR-FIELD RADIATION PATTERN 


Poe purpose of this far-field radiation pattern, the helical antenna 
axiS 1S equivalent to the Ex = 0 grid line. The pattern is essentially 
symetric when rotated about the antenna's axis. 


Ex(@) := fo. 4 cos( 0) Ey(0) := fo 2} si 


0.98 
0.76 


0.54 


Za 


Ey(®) 0.1 Roky Za 
=O) 2 
0.34 


—0.56 


—0.78 





—0.8 0.64 —0.48 0.32 -0.16 5.551 10 


Ex( 6) 
nie hile 


a7 


PARABOLIC REFLECTOR ANALYSIS WITH HELICAL ANTENNA FEED 


The helical antenna Mathcad application (Reflector Opticon) will eeenemee 
the following additional parameters (Note: a 2 00re22 wil eto Mtouwemes 

secondary field parameters mime rns application, and inorder to speed up 
computation time Corresponding vectors are created with subscript eigen 


Px,y = E-Field Integral and Summation Equation Vector Components 
Eax,ay = Aperture E-Field Components 

0',p = E-Field Integral and Summation Equation Parameters 
Eg = Secondary Field Theta E-Field Cempeowene 
E4 = Secondary Field Phi E-Field Component 

U = Secondary Field Radiation Intensity 

Praq = secondary Field Radivated Power 

D, = secondary Fields tirecremyvsse,, 

G = Secondary Field Gain 

fap = Aperture Efficiency 

fé- | = “Spillovers tirerency 

fe «= @ laper Etre eney 

i = Phase Efficiency 

& = Polarization Erticienay 

Eohmic = Ohmic Efficiency 

Ag >= Retlector Eitecervegnren 

he = Retlector Ef tecervesieitgnr 

EIRP = Reflector Effective Radiated Power 
Ex,y,z = Secondary Field Cartesian Components 
6, = Reflector Unit’ Polarizagmeweiecctor 

PLF = Reflector Polarization Loss Factor 

Imin = Minimum Dirstanceseo Secondasyeearywela 
0, = Reflector Periphery Coaltitude 

Or bp = Theta and Phi for Polarization Loss 


The following additional data must be inpueEted based on known or 
estimated data: 


Eo -> Electric YFieldsscale, Factor 

Q@ref = Radius of the Mouth of the Reflector 

foc = Reflector Focal Length 

& = Blockage Efficiency 

Esp = Spar Efficiency 

tl,t2 = Radiated Power Increments for Phi and Theta 
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il = Secondary Field Increments 


N,M = Summation Increments 

oes >econdary Field Observation Distance 

X,Y,Z = Polarization Loss Factor Coordinates 
Weeeameuacoming Wave Unit Polarization Vector 


Enter additional input data here: 


foce=> (meters) aref := 5 (meters) 


(Note: (£) should be greater 
than minimum distance to the 
far-field for the feed antenna) 


eb = .96 (dimensionless) 

esp = .95 (dimensionless) tl =10 (increments) 
t22= 10 (increments) i] =10 (increments) 
N =10 (increments) Mir=al0 (increments) 
ff2 = 10° (meters) x = 1000 (m) 

J _ 

ao y = 1000 (m) 

f2 
ow2 =] | (dimensionless) 

fe z =1000 (m) 

0 


ies, 


Establish@integrammand Summation Increments: 


(radians) 


(radians) 





92 =-(4]- 10°, Se ea eer Oe 
36 11-18 36 36 


(m) 


(m) 


Calculate reflector geomet rie Patameters: 











Qo i= 2 aan sy (radians) 
2-foc 
60 = 0.9273 (radians) 
4-foc” + rp 
p( rp) := (m) 
uy 4. fi 
ae rp | 
O@p(rp) ‘= 2 tan | (radians) 
2-foc 
4. foc? + (rp!) 
pl = (m) 
4-foc 
6p :=2-at ae (radi ) 
= 2-atan radians 
= ot 





1 
= 
0 36 
Set lees 
r:=10-aref 
_ ., S-arefy 
n 


rmin ‘= max(r) 


rmin = 2.0014 10° 
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(increments) 


(increments) 


(radians) 


(radians) 


(radians) 


(radians) 


(m) 


(m) 


(m) 


(m) 


(m) 


@eeoeure Electric Fields Baxs Bay! 


Eax(rp, $l) := (E@( 8p(1p))-- cos( $1) + Ed(@p(rp))-sin(g1 ))-e) SC PUP)= foe2) 
Eay(rp,61) = (EQ 6p(rp))-sin($1 ) + Ed( 6p(rp)):cos(1)):-e! S¢PCP)= foe?) 
Eaxl , ‘=| B0(@p!,) -(- cos(11,,}} + E6(@p1_,} sin($11,) |. Kk P1q- foc} 
Eayl, ,:= (E0(@p1,) -sin(411,) + E¢(Op1,) cos(911,))-¢ (Pla foe} 


Pegeertc Field Integral Equation Vector Components Px, yi 











aref 
Ax(62,62,61) =| Eax(rp, 1 )-e -k-rp-sin( 2 )-( cos( $1 )-cos( 2) -+ sim $1)-810( $2)) ey 
0 
2:7 
Px( 62,62) :| Ax( 82, 62,61 ) dol 
@) 
aref | | | | 
Ay( 82,62 ,61 ) :| Eay(rp, 61 )-e! Bsc C2) (rere (PU) coat #2) esi QE) S82) dep 
0 
2:1 
2.42) | Ay( 92, 62,91 ) dol 
0 
(223) -(onlrig)-con(o29) sing) sin( 229) 
j -k-rpl-sin(622 a) -(cosl$11p) -cos(622p) + sin{o11p) -sin{622 
Axl, 5 = Eaxl, ,¢ ‘rp| drp 
‘0 
2:1 
Pxl ap = Ax], , 0611 
“0 
oie (23) (corn) connag) + sinorig) sin(n29) 
j -k-rpl-sin(€22 4] - {cos {11p) -coso22p) + sin{o11p) -sinlo22 
Ayl, 3° Eayl, rp! drp 
0 
2-1 
Be Ayl, , 4611 
!0 


red 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


Secondary Field Electric shield Components Eg, 


oe -j -k-rff2 
592(02,62) 22 K2 (px a2, 42)-cos(62) + Py(02,62)-sin(62)) 
-n-rff2 
faker one 
E62(62 ,62) -= -(Px( 62, 62)-- sin(o2) + Py(@2, $2 )-cos( $2 ))-cos( 62) 
on 
jo emaee | 
£03, 2 re (Px1, , 608(622, + Pyl, ,-sin(622,)] 
os 
i eae | 
E93, 2p = — [Pall es (- sin($22,) + Py], ,cos( $22, } cos(022,) 
= 


Secondary Field Radiation Intensircy jO fe 


Pe 
U2(02,42) — 1592(62,62)1)7 + (1E42(02,62)1 7] 
. 
rff2? 
U3,5 °F [ (LE, al) *+ (EPs) 
0 





Secondary Field Radiated Power Prad: 


Prad2 := (2 (4) 2 20° aob22,) sin (022 4 


Prad2 =0.19318 


Secondary Field Directivit Da: 





_ 4... max( U3) 
Prad2 


Do22: 


Do22 = 1.04762 10° 


Eo 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(W/solid ang) 


(W/solid ang) 


(W) 


(W) 


(dimensionless) 


(dimensionless) 


Seeeeover Efficiency es: 


Bo 
| U(@p):sin(@p ) dp 
LU is 
t 
U(0p )-sin(@p ) dp 
“0 
es = 0.19013 


Peteeeeetficiency ef: 


Go 2 
2 
eet an 2p dép| | -2-cot] |— 
Prad 2 
9 oll 
0) 
SOP) AT sin( Op) dop 
Prad 
“0 
et = 0.50837 
Pease Efficiency Epi 
-4.-6:10 >-aref 
re A 
ep = © 
ep = 0.68574 
Beerture Efficiency fap. 
ex := 98 (dimensionless) eohmic ‘= 98 
gap ‘= eS t-ep-ex eb: esp-eohmic (dimensionless) cap =0.05805 
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(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Effective PSOErOpIC Radiated Power EIRP: 


EIRP22 := Prad2:Do22 (W) 

EIRP22 = 2.02379 10° (W) 

Gain G: 

G22 :=eap:-Do22 (dimensionless) G22dB ‘= 10-log(G22) (dimensionless) 
G22 =6.08192 10° (dB) G22dB = 37.84041 (dB) 


Effective Aperture AL. 


Ae22 := eap-n: aref” (mé ) 


Ae2?2 = 4.55961 (m7) 


Radiation Resistance Rye 








Rr22 =2. a (Q) 
([Io] ) 
Rr22 =0.38636 (Q) 
Effective Heighe Rae 
he22 ‘=2. IRr22 Ae42 (m) 
No 
he22 =0.13672 (m) 
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Polarization Loss Factor PLF: 


nix? + : (radians) ! 
6p2 := atan y tp2:=atan(2 (radians) 


Zz 


6p2 = 0.95532 (radians) dp2 =0.7854 (radians) 
E64 := E€2(0p2, ¢p2) (V/m) E94 := Eo2(8p2, op2) (V/m) 
Ex2 ‘= E64-cos( 6p2)-cos( $p2) — Ed4-sin(¢p2) (V/m) 

Ey2 := E64-cos(6p2)-sin(¢p2) - E¢4-cos( op2) (V/m) 


Ez2 := - E@4-sin( 8p2) (V/m) 


Ex2 


_  _ ae Ey2 (dimensionless) 


J [Ex2| )? + (|Ey2| )?+ (IEz2|)? 4Ez2 


0.12034 + 0.48533) 
oa2 =| 0.12034 + 0.48533) (dimensionless) 
0.60556 - 0.36501) 


——)\ 2 
Rer2 := el) (dimensionless) 


PLF2 = 0.25003 (dimensionless) 
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PARABOLIC REFLECTOR WITH HELICAL FEED FAR-FIELD ELEVATION PATTERNS 


E-FIELD RADIATION PATTERN: 





Q2 ‘=. <|- i0° 2h 2 Go eee 
i118 36 36 


ole + ;| 
2 


X2(92) ‘= 














E@2 ( 3 
z 


0.032 


0028 


0 024 


0.02] 


¥2( 82) 0017 


0.014 


0.0] 


0.003 





0 
—0.003 —0.002 =0.002 —0 001 “0 001 0 0001 0001 0.002 0.002 0.003 


— “PREY 
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H-FIELD RADIATION PATTERN: 


P02 ) -= [Eo2( 02.0 )| cos} 2 + kj Y3102) := |Eo2(02.0 )| Ws + a 
a a 


0.035 


| 0.032 | 

| 

( LE 
v 024 


0021 
} Y3( 62) 0.017 
0.014 


0.01 





0 
—0.003 —U 002 “0002 —0001 —0 001 0 0.00) 0.00} 0 002 0.002 0.005 


ee eel 
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THE SPIRAL ANTENNA (RE? LECT Or SOR eaee 
MATHCAD SOFTWARE-DISH SPI. MeD 


Spiral antennas are a family of two or three dimensional devices that 
possess frequency independant parameters over a wide bandwidth. Spiral 
antennas are commonly used for direction finding, satellite tracking 
and missile guidance. 


The planar spiral may be of the Archemedean, log-spiral, orvequian la. 
type. All three radiate two main, circularly polarizZeg@lioves 
perpendicular to the plane of the antenna. Additional gain for planar 
spirals may be achieved by placing a metal cavity on the side of the 
antenna with the unwanted lobe. The cavity may be empty or be filled 
with electromagnetic energy absorbing material. |these appitesreieane 
principly examine the equiangular planar spiral and do not account for 
cavity backed effects. 


(Note: Mathcad equations cannot use symbolic subscripts. Therefore, 
symbols like »% will immediately follow the parameter in equations. ) 


The spiral antenna Mathcad applicatwens will compute the following 
pabamieters fOr equiangular plana Spatals: 


k = Wavenumber 

i = Wavelength 

Do = Direectrvicty 

Fo = Electric Field Component 

U.~ = /.Radiaticn Imeem re, 

Umax = Maximum Radiation Intensity 
Prag — Radtaced rower 

G = Gain 

EIRP = Effective Isotropic Radiated Power 
Aen = Maximum Effective Fperture 
BW = Bandwidth 

Emin — Minimum Distance (our 4u amelie 


Ry = Radiation Resistance 


Nem = Maximum Etfective sievang 


fhigh = Upper Operating er zequeme, 

fiow = Lower Operating Frequency 

a — AMy port, onmene nth edge of a Spear 

Eex = Equiangular Planar Spiral ExpansiompRat io 
Z1 = Planar Spiral Input Impedance 

[ = Voltage Reflection Coefficient 

fry = Ketlecticon Eriievency, 


os 








fee POL aALIZaALLOCON Loss ractcy 
high = Upper SCperacing wavelencth 
Bem UOWer Operating Wavelencth 


Peeaeeeeeranar Spiral Electric Fieid Amplitude 


The followino data must be input based on known or eStimated data: 


M = Mode 

N = Number of Spiral Arms 

Meomocecuiency Of Interest 

1 = Number of Increments for Far Field Radiation Patterns 
iggeeeeiastance of Far-Field Calculations 

Ty 


Seinput Current at Antenna Terminals 
ieee raral Feed Point 
a = Flare Rate 
$54] = Angular Arm Width of nt) spiral Arm 
®, = Azimuth to Compute Expansion Ratio 


Peeomeeoverall Radius 
Ppeeeeoource Strength Constant for Planar Spirals 


Gemeeewave Unit Polarization Vector 
Memeecnaracteristic Impedance of Feed Assembly 


eee agutengular Planar Spiral Unit Polarization Vector 


THE PLANAR SPIRAL ANTENNAS 


Enter input data nee 


Le (arms) f-=3.5.10'° (HZ) 
M ‘= | (mode) 

lo :=] (amps) 
mioce: Mio. is N-1) 
f= 50 (ancrements) rff-21 10° (meters) 
a= 2 | (dimensionless) ro = OO! (m) 
Seed fan Eo := 10° (V/m) 


189 


| 
f 
e 
6:= (radians) ow =] } (dimensionless) 
. = 
J 
eset ‘“ 
5 poy 


Lo = 100 (Q)) 


Galculate pianar spiral antenna Geometric Parameters and define 


constants: 


¢ :=29979 10° (meters / sec) Ny = 120-4 Bel 
c 
‘S— 2 i 
/ : (meters / cycle) 6, | or? (Parads 7am 
SOF 
» i 4s. yee 
2 =856543°10 > (meters / cycle) Wy za lo! Cea 
] 
Zz 
(dimensionless) 
oa:=] | 
i 
Oil) 
Calculate planar spiral antennas corameue wom 
Def inewancu lar ost oecmem 
EGOMGy seman lon 
l 
Se) Oe ee >- 10° = (radians) j= Oe ean (radians) 
ies ] 


oo Cod ed 


io 'O 


Bueoomce tO Far-Fiele c,.-: 








ies | Os (m) nen (m) 
sR? 

IT, — (m) 

mK 
rin := max(rr) (m) 
rmin =0.] (m) 
Wavenumber kK: Radial Distance to nth Spiral Edge ©: 

2 Or— 6 

kia (m7 1) Pe core i ) (m) 
} 

9733552 10° tm 2) 1.22) =4 0091710 > (m) 

Expansion Ratio fy: 
eex(n, or) Secu + 2) (dimensionless) 
r(n, or) 
eex( 1,2 2) = 4.00917 (dimensionless) 
Bandwidth BW: 
Equiangular Spiral EOgaeer1eodic Soir al 
Ahigh ‘= 4-ro (m) Ahigh] := 20-ro (rm) 
Cc a 
fhigh := fhigh]) :=——— 

7 Ahigh aac ° zhigh| (Hz) 


ial! 


10 


thigh =7 49475 10 ee thigh =~ 49475 1" he 





dlow :=4 R (m) 
flow ‘= : (Hz) 
low 
flow = 7 49475 10° (HE) 
BW .= fhigh - flow (Hz) BW] .= thigh] = flow (Hz) 
BW =6 74528 10!° (Hz) BW] =749475+10° ie 
Electric Field E(0$) and Electric Brel teude een 
wieQ0 1 (increments) 
; ore p\M 3 ss (v4 2) — na 
Eo-k*-cos(6) (1+) a cos(@)) : on (2 € y 
=) 
E(6@) ‘= : — (V/m) 
sin(@)” rf¥ 
9 M (* atan( a cos(9)) 
cos() tn et 
9) 
A(8) rr ae (V/m) 
sin(@) ar near cos( 6)” 
M 
w( 8) := a int +a" eon ai) + atan(a-cos(6) ) (dimensionless) 


pagar e| 


nx 
M (-+2}-1] (V/m) 





J . 


E48) = (8)? 4 = 
rf 


1g 2 


n Ww 
ee oe (SM) tan{a cos(- 3 4 =) 
n Ww Z 1 a 2 i 
col 3 | {an|—__——_J] _-e (V/m) 
Aan i= 2 1 2 
Ww 
2 
ip ee Ww 2 n Ww 
sinj-— + m— | [1+ a™-cos|-— + — 
| 2 : | : 
Peag@edtron Intensity U(8): 
= .(acoy’) 
U(6) :=——- \A(6) (W / solid angle) 
2-1 
Ul =—|{a1,}] (W / solid angle) 
Ww 2, 
Umax ‘= max( U1) (W / solid angle) 
Umax = 2.2285°10 > (W / solid angle) 


Radiated Power Prad: 


2 

Prad ‘= 4-2 U(6)-sin( 8) dé (W) 
0 

Prad = 7.94528°10- (Ww) 


Emecectivity Do: 


0 - 4x Umax (dimensionless) 
Prad 
D. = 3.52462 (dimensionless) 


eo 3 





Radliaticr Resistance Rr: LASUS SimMseecarnce 2a 
2 Prad | 
=e (Q) ,. N30a 
2 At ——$——_— (02) 
([Io] ) M 
Sinn. 
N 
Rr = 0.01589 (2) Zi = 1 88496" 10° (Q) 


Voltage Ret lect vomiccen hialememis Reflecton EEElecilency ets 





reir (dimensionless) om :=1- (rl) (dimensionless) 
Zit oo 
l = 0 30675 (dimensionless) ay =0 9059] (dimensionless) 
Galn G: 
oie tans De (dimensionless) GdB ‘= 10 log(G) (dB) 
G =3 19298 (dimensionless) GdB =5 04196 (dB) 


Effective TPSOtropic Radiated Power (EIRP): 


EIRP ._ Prad G CW) 
= 
EIRP =0.028 (W) 


Polarization LOSS Factor (2Lr i 


S\-2 
PLE = (low cal) (dimensionless) 


PLF = ] (dimensionless) 
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j 


Mere rmerm ES fective Aperture & 


He D, av PLE 
Aem (= ———————-——$ _—— 
4-4 


Aem = 1 86416 10> 


Memmum Effective Height <h,,) : 


Rr Aem , 
No 


hem °:= 


hem = 5 6063°10 > 


oS 


THE EQUIANGULAR PLANAR SPIRAL ANTENNA FAR-FIELD ELEVATION PATTERN 


For the purpose of this far-field radiation pattern, the sorrel ano 
lies parallel to the Ey = 0 grid line and 1s centered at the origin. 
The magitude of the electric field pattern 1s rotationally symmetric 
With respect to the Ex=0 grid line. The equiangular planar Spiral 
antenna possesseS a Mlrror image radiaton pattérn in the -yY Maweee 


E$x(6) = | Eoc 6) | cn 4 Egv(G) := |E4(6)| snl 4 


0.001 

ACC 
an || 

Eoy( @) tt tLe ELLE 


Svea 





v= FIELD 
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PARABOLIC REFLECTOR ANALYSIS WITH SPIRAL ANTENNA FEED 


ieee prtal antenna Mathcad application (Reflector Option) will compute 
Meemeeol Owing additional parameters (Note: a 2 or 22 will follow most 
Seeemcdary field parameters in this application, and inorder to speed up 
pem@emedtiOn time cOrresponding vectors are created with subscript A & B: 


Px, y pmeeeriecla Integral and Summation Equation Vector Components 
Bax,ay =weepekrteure E-rFiela Components 

Cpe — ) ©-Fielda Integral and Summation Equation Parameters 
Pome ccondary Field Theta E-Pield Component 

E 4 =e secondary Field Phi E-Field Component 
Meo secondary Field Radiation Intensity 

Peep Secondary Field Radiated Power 

See oecondary Field Directivity 

G = Secondary Field Gain 

fap = Peerture BrLficiency 

Gees oPillover Efficiency 

Mame laper Efficiency 

oo) ir Phase Efficiency 

WumeeeeerOlarization Efficiency 

ates | |OOnmic Efficiency 

Demeeeretlector Effective Area 

Mupeceenet lector Effective Height 

EIRP = Reflector Effective Radiated Power 

Ex, y,Z = Secondary Field Cartesian Components 
amet lector Unit Polarization Vector 

Peres Reflector Polarization Loss Factor 

See «Minimum Distance to Secondary Far-Field 
[eee tlector Periphery Coaltitude 

O51 bp = Theta and Phi for Polarization Loss 


The following additional data must be apueted based on known or 
estimated data: 


Sof Radius of the Mouth of the Reflector 

foc = Reflector Focal Length 

fe SloecKage Efficiency 

Esp = Spar ELilevency 

t1,t2 = Radiated Power Increments for Phi and Theta 


ey 


1l = Secondary Field Increments 


N,M = Summation Increments 

Yes> = Secondary FielaeObservartion is tance 

X,Y,Z = Polarization Loss Factor Coordinates 
Oo, = %Incoming Wave Unit, Polar zations Veccra. 


Enter additional input data here: 


aref ‘= 0.32 (meters) 
foc ‘= 0.147 (meters) 
esp := .95 (dimensionless) 
2 := 10 (1mcrements) 
N ‘= 10 (increments) 
rff2 := 10° (meters) 
ie 
ya 
OW2 := (dimensionless) 


eb := 96 
tl -=10 
1] := 10 
M := 10 
x ?= 1000 
y := 1000 
z >= 1000 


(dimensionless) 


(ancrements) 


(increments) 


(increments) 


(m) 


(m) 


Establish tncegral and Summation Increments: 


(radians) 


oie 0 —..2-n 
t] 


92 2103 274 103. 2-24 10° 
1] 





Pees io —— - — - 10°. = - 10° 
36 i] 18 


36 36 
+0, aref (m) 
t 
ipl,= aref-— (m) 


(radians) 


A:=1.N 
B:=].M 
B 
i= 2 — 
: M 
2 =i 
M 


Calculate reflector Jeemeu@rer parameters: 





Qo := aan | (radians) 
2:foc 
@o = 1.65544 (radians) 
4-foc? + rp’ 
orp) = om 
-foc 


(radians) 


@p(rp) ‘= 2 aan| | 


‘foc 


: 4-foc” + (rp!,)° 





pl: no 
= 4-foc 
rp 
pl, <2 3 (radians) 
2: foc 





=o? 
Gor 10. aref 
_ ., S.aref 
é d 


rmin ‘= max(r) 


rmin = 95.64028 


oo 


(increments) 


(increments) 


(radians) 


(radians) 


(radians) 


(radians) 


(m) 


(m) 


(m) 


Aperture Electric Fields Eames Eays 


Eax(rp, 6] ) ‘= EX @p(rp))-sin(g] )e) “(PUP = fe?) 


Eay(rp, $1) “= Ex @p(rp))-cos( gl )-e2 S (PUP) = for?) 


Eaxl , ,, = E¢(6p!,,)-sin(¢l Ip)e ‘k: (pl = foc 2) 


pl — foc-2) 


Fayl, 's E4(6p1 ,}-cos(¢1 In)e b( 


Electric Field imeeg val Equation Vector Components Px, Ye 


aref 


ance. 24n | Eax(rp, $1 )-e k-rp-sint 2 )-( cos( ES SAN CHIC) 
0 


Uy, 
pu. =| Ax( 62, 62,61 ) do! 
0 


aref 
Ay(62, 62,1) :| Eay(rp,o1)-e) EP sin( 82) ( cost 41 )-cox $2) 4 sin ¢1)-sim 02), 
0 
20 


p22) =| Ay(82, 62,41 ) dol 
0 


Fax] ae k-rp]-sin(022,4) (cos(41 1p) -cos($22p) +sin((1tp)-sin(@225)) 


200 


drp | 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


aa ce al 


— 


Secondary Field Electric Field Components Eg, Pot 


ey kif? 
E62( 62, $2) Ee EE px 62) 62). c0s( 42) 4 Pv( 82 , $2 )-sin( $2 )) 
2-n-rff2 


ed kar 
E62( 62, 42) = 4 SS (Px 02,42) sin 42) 4 Py( 02 , $2 )-cos( 2) )-cos( 02) 


1: & 


E03, , =. (Px1, , cos(§22,} + Pyi, ,-sin(622,}} 


E¢3. 5 eee Pct,» (- sin(422,}) + Py] a.p005(422,) | cos( 622, } 


Secondary Field Radiation ines ey \o2 


2 
U2 62,42) = {k |£02(02,42)1)°+ (1E4202.62)1)7] 
” 


exif2” 
U3, 3 = | ({ E83, .|)°4 (E43, |)’ | 





Secondary Field Radiated Power Prad: 
2-0 nN 
— —_— | ee « 2g) wm . 2 
Prad2 | } - | ; (2 U2(022,.42 a) sin(9 ) 


Prad2 = 1.84163 10 |! 





Secondary Field Directivit Dee 


Do?? := 4-nmax( U3) 
Prad2 


Do22 = 4.59318 10° 
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(V/m) 


(V/m) 


(V/m) 


(V/m) 


(W/solid ang) 


(W/solid ang) 


(W) 


(dimensionless) 


(dimensionless) 


spillover Elrtevenevmcc 


Go 
| U(@p )-sin( 6p ) dép 


| ig? 
gs .= 
se Os 
2 
U( 8p ): sin( 8p ) d6p 
aCay 
es = 1.00042 


Taper Efficiency ee: 


2 


2 
60) 42 an d@p 2a (&)] 
; Prad 2 Z 


Bo 


et = 0 45135 


Phase vb Et rcteney bp 


ep = 0.97222 


Aberture ,ELtVcrency fap. 


ex -= .98 (dimensionless) 


cap -= eS et -ep-ex eb esp cohmic (dimensionless) 
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eohmic .= 98 


cap = 0.38452 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


(dimensionless) 


Brerective Heoenoere Radltated Power EIRP: 


Elk? 22 := Prad2-Do22 


EIRP22 =8.45893°10 ° 


Gain G: 


G22 ‘= eap Do22 (dimensionless) 


G22 =176615°10° (dB) 


Effective Aperture ea 


Ae22 *= cap 1: aref- 


Ae22 =0.1237 


Radiation Resistance Ry: 


Prad2 
([Tol )? 


Rr22 :=2 





Rr22 = 3.68326 10 |! 


Effective Berghe Ne: 


Ae22 
No 





he22 :=2- [Rr22. 


he22 =2.19868 10 ” 
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G22dB ‘= 10-log( G22) 


G22dB = 52 47028 


(W) 


(W) 


(dimensionless) 


(aB) 


(Q) 


(m) 


(m) 


Polarization LOSS Factor PLE: 


7 x’ + y’ (radians) = yy 

6p2 := atan| ——_—__—_ op2 := atan|— (radians) 
= 
6p2 = 0.95532 (radians) ¢p2 = 0.7854 (radians) 
E64 := E62(6p2, 6p2) (V/m) E64 := Eo2(@p2, ¢p2) (V/m) 
Ex2 := E64-cos( @p2)-cos( ¢p2) — E¢4-sin( $p2) (V/m) 
Ey2 ‘= E64-cos(@p2)-sin($p2) - E¢4-cos($p2) (V/m) 
— Sj 
Ez2 ‘= - E64-sin( 6p2) (V/m) 
x2 

oad = Ey? (dimensionless) 





({Ex2})?4 C{Ey2l 24 C{Ez2I)? \Ez2 


0.40672 + 0.03531) 
oa2 = |-0.40672 + 0.0353]; (dimensionless) 
0.81344 - 0.07063) 


PLF2:: eile (dimensionless) 


PLF2 =0.16667 (dimensionless) 
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PARABOLIC REFLECTOR WITH SPIRAL FEED ANTENNA 


E-PANE RADIATION PATTERN: 














il :=20 os [2)- Mig eres Oe 
36 1] 18 36 36 

X1(02) ‘= pin(an,4| cos(tr+ 2 Y1(02) ‘= e02(o2.2) sin(o2+ 3 

2 2 2 2 
3.09 10 /° 
3.08910 !° 
3.08710 !° 
3.086 10 !9 
3.08a¢10 19 
-15 


Y1( 82) 3.083 10 
g !5 


3.082° 1 


3.08 10 !9 


3.07910 !° 


3.077 10 


3.076 10. 





a PED 
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H-PLANE RADIATION PATTERN: 


X2( 02) ‘= (£62(62,0)| cos 7 ;) Y2(@2) := (£42( 62 ,0)| oof 2 


“ |) | fete ates 
7 Eat 1] ah t 
1.06910 © al | 


i Sea RLEEE 
8.2 (ume 


Yx@2) 7°10. 


; cote 

LANGA | 
SERRARP AREER 
lel scleelselale lal la 


=6)10-¢8 108-36 108-2 4108-12 ue 0 1210 ~ 2.410 ~ 3610 ~ 48 10 6 10 © 
X2( 82) 


~J 

— 

Laienaanll Ey 
—— 


a 


34100 


uf 


= Je hIELD 
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EO. 


oA. 


a2. 


7. 
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